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PREFACE 


This  technical  report  covers  the  work  performed  by  NESCC 
during  the  period  May  1,  19bii  to  March  31,  1963,  pertaining  to 
the  determination  of  damage  In  the  lower  Chesapeake  Bay  due  to 
long  waves  generated  by  underwater  nuclear  explosion. 

It  is  emphasi74ed  that  this  work  has  required  the  development 
of  new  theories  relevant  to  the  hydrodynamics  of  long  waves. 

These  theories  may  find  many  other  applications  than  those  which 
have  been  carried  out  for  the  particular  problem  under  study.  On 
the  other  hand,  because  of  the  timing  of  the  project,  some  of  these 
new  theoretical  tools  have  not  been  fully  exploited  as  they  could 
have  been  for  obtaining  more  accurate  and  reliable  quantitative 
results  in  the  lower  Cheaapeakc  Bay  areas. 

Despite  these  new  theoretical  developments,  and  due  to  the 
complexity  of  the  problem  under  study,  a  number  of  gaps  remain. 
Hence  some  practical  results  rely  more  on  engineering  judgment 
and  approximate  methods  than  on  scientific  calculations.  New 
theoretical  or  experimental  research  will  certainly  permit  a  better 
quantitative  investigation.  However,  it  Is  believed  that  the  general 
conclusions  presented  in  this  report  will  not  he  greatly  modified 
by  increasing  the  scope  of  this  study,  but  they  can  be  given  with 
more  accuracy  and  details. 

The  work  on  this  project  was  carried  out  in  the  NESCO  office 
at  Pasadena,  California,  Washington,  D.  C. ,  and  Houston,  Texas. 
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The  work  in  the  Pasadena  office  was  under  the  supervision  of  Dr. 
Basil  W.  Wilson,  with  added  contributions  from  Dr.  James  A. 
Iir-?>drickson  and  Miss  Lois  Webb.  The  work  in  the  Houston  office 
was  under  the  supervision  of  Dr.  John  C.  Freeman,  Jr,  ,  with 
add  contrihut.ion.s  from  Dr.  Larry  Armijo  and  Mis.?  Mary  Ann 
Noser.  The  work  in  the  Washington  office  was  under  the  super¬ 
vision  of  Dr.  Bernard  Le  Mehaut^,  who  was  also  project  principle 
investigator  and  coordinator,  with  added  contributions  from  Dr. 

J.  Ian  Collins,  Mr.  Raniinder  S.  Grev/al  and  Mr.  Carl  Fast. 

Dr.  Le  Mehaute  is  credited  with  the  content  of  Volume  I 
of  the  report  dealing  with  the  engineering  aspects  of  the  study. 
Various  members  of  the  NESCO  team  have  made  original  contri¬ 
butions  to  the  advancement  of  hydrodynamic  theory.  Their  scien¬ 
tific  contributions  are  reflected  in  the  various  appendices  of  this 
report  presented  in  Volume  II: 

Appendix  I:  Long  Waves  Generated  by  Nu  .  lear  Explosions; 

Type  of  Waves  and  Initial  Decay  in  Deep  Water, 
by  Dr.  Basil  W.  Wilson. 

Appendix  II:  Theoretical  Considerations  and  Computations  for 

Water  Waves  Produced  by  Explosion,  by  Dr. 
T.arry  Armijo  and  Mis.s  Mary  Ann  Noser. 

Appendix  III:  Surface  Waves  Generated  by  a  Disturbance  on 

Sea  Bed  in  Constant  Depth  Open  Sea,  by  Dr. 
James  A.  Hendrickson. 
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Appendix  IV ■.  'I'ho  F^vinciple  of  Supe rpoBitinn  and  fhe  Theory 

of  Caxir.hy -  PoisKon  by  Dr.  Bernard  Le  Mehaut6.. 

Appendix  V:  The  Shoaling,  Diimping,  Breaking  and  Run-up  of 

Long  Waves  Over  the  Continental  Shelf;  On 
Saturated  and  ‘Nonsaturated  Breakers,  by  Dr. 
Bernard  Le 

Appendix  VI:  The  Wave  Riin-upby  the  Method  of  Characteris¬ 

tics,  by  Dr.  Bernard  Le  Mehautc.  Dr.  John 
Freeman  was  scientific  advisor  and  Mr. 

Ramindcr  Grewal  did  the  immerical  computations. 

Appendix  VII:  Two-Dimensional  Nonlinear  Wave  Motion  in  an 

Estuary,  by  Drs.  John  C.  Freen^an,  Jr.  ,  and 
Larry  Armijo,  and  Miss  Mary  Ann  Noser. 

The  undersigned  would  like  to  express  his  appreciation  for 
the  assistance  rendered  to  the  project  by  Dr.  Lewis  McCammon, 
Associate  Director  of  Engineering  at  the  Pasadena  office  of  NESCO, 
and  to  Mr.  ZemgaLs  for  supervising  the  drafting  of  figures. 

Acknowledgement  is  given  to  Mr.  Kenneth  Kaplan,  of  United 
Research  Service,  Inc.  of  Burlingame,  California,  who  furnished 
the  input  data  for  the  project  and  acted  as  a  consultant  to  NESCO. 

Charles  L.  Bretschneider 
Director,  Washington  Office 
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ABSTRACT 

This  report  deals  with  the  effect  of  tsunami-type  gravity 
waves  caused  by  huge  underwater  disturbance,  on  the  shoreline 
environment. 

In  particular,  the  cases  of  three  different  power  disturbances 
off  the  United  States  Atlantic  Coast  and  their  effects  in  the  lower 
Chesapeake  Bay  are  analyzed.  It  is  shown  that  because  of  the 
natural  protective  action  of  the  continental  shelf,  a  maximum 
possible  amount  of  damage  is  reached  for  a  given  power.  If  a 
disturbance  is  due  to  a  larger  power,  the  difference  of  energy  is 
dissipated  over  the  continental  shelf  by  spilling  breakers. 

The  flooding  areas,  dynamic  destruction  and  harbor  reson¬ 
ance  due  to  this  optimum  power  disturbance  are  calculated  and/or 
evaluated. 


xxi 
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PART  A 


PURPOSE  OF  THE  STUDY 
AND 

GENERAL  CONCLUSIONS 


by 

R.  Le  Mehaut^ 
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CHAPTER  I 

PURPOSE  OF  THE  STUDY 

From  time  to  time  disturbances  at  sea  produce  long  waves 
which  cause  considerable  damage  to  populated  areas.  These  waves, 
so-called  tsunami  waves,  are  most  often  due  to  a  disturbance  of 
seismic  origin.  The  time  has  now  arrived  when  such  destructive 
waves  can  also  be  generated  by  explosions  of  large  nuclear  devices. 
Therefore,  in  view  of  the  range  of  wave  heights  to  be  expected  from 
such  explosions,  impetus  has  been  added  to  the  immediate  need  to 
study  this  problem. 

The  first  purpose  of  this  study  is  to  determine  whether  or 
iioi,  ilieie  is  a  prouabliily  oj  daiiiage  along  coastal  areas  due  to 
large  yield  nuclear  explosions  at  sea.  Then,  if  this  probability- 
does  exist  (and  it  has  been  found  that  it  does  exist),  the  study 
would  include  an  estimation  of  the  damage  to  coastal  and  harbor 
facilities  due  to  these  tong  waves  in  tlie  lower  Chesapeake  Bay,  and 
possible  engineering  solutions. 

The  size  of  the  assumed  initial  disturbance  at  the  free  surface 
and  the  wave  trains  due  to  an  explosion  for  a  point  twenty  miles  from 
the  disturbance  were  provided  to  NESCO  by  United  Research  Ser¬ 
vices  for  three  cases  corresponding  to  three  different  bomb  yields. 
For  purposes  of  this  report  the  smallest  yield  is  defined  as  Case  I 
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and  the  largest  as  Case  III  (see  Fig.  1).  The  origin  of  the  dis¬ 
turbance  is  assumed  to  be  some  370  nautical  miles  from  the 
Atlantic  coast  (Fig.  Z),  Tlie  evolution  of  these  wave  trains  is 
then  calculated  from  ZO  miles  from  the  explosion  to  the  extreme 
end  of  the  estuaries  and  harbor  basins  until  their  effects  become 
negligible. 

A  practical  analysis  of  the  travel  of  the  wave  trains  in¬ 
volves  a  succession  of  various  hydrodynamic  problems  to  be 
solved  by  completely  differing  methods. 

The  first  problem  consists  of  the  calculation  of  the  evolu¬ 
tion  of  the  wave  train  in  deep  water  up  to  the  toe  of  the  continental 
slope.  Then  shoaling  and  reflection  effects  due  to  the  continental 
slope  have  to  be  evaluated. 

tr a  shelf  are  snhiecteH  to 

shaoling,  refraction  and  relative  damping  by  bottom  friction.  Ul¬ 
timately  they  break  up  into  a  succession  of  spilling  breakers,  then 
bores,  and  finally  run-up  over  the  coastal  dunes  causing  flooding 
and  destruction  along  great  stretches  of  the  Atlantic  coast.  Also, 
when  the  waves  penetrate  river  estuaries  and  bays,  they  diffract 
and  give  rise  to  multiple  reflections  causing  severe  local  agitation 
combined  with  wave  run-up  and  flooding.  Finally,  upon  reaching 
harbor  basins  they  are  the  cattse  of  great  agitation  induced  by 

r.  -i-.  i  1  1  1  f  #  -  /-•  K  m  •»*  o  r*  I  o  V  i  CJ  I  i  » •  C  ea  Q  c> 
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Qualitative  details  of  all  these  different  aspects  are  described 
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FIGURE  I 

ASSUMED  LOCATION  FOR  THE  EXPLOSIONS 
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in  Part  B.  Each  of  them  has  been  sufcjjected  to  theoretical  studies 
and  applications .  Great  use  has  been  made  of  existing  theoretical 
and  experimental  studies.  However,  it  was  found  almost  impos¬ 
sible  to  provide  a  satisfactory  answer  to  these  complex  problems 
without  an  extension  of  existing  knowledge.  Hence,  now  research 
on  these  problems  bccarnc  a  major  part  of  this  study. 

Some  of  these  new  theories  can  be  further  refined  and  im¬ 
proved.  Others  which  will  be  directly  related  to  the  problem  under 
study  have  not  been  investigated  in  this  report  because  of  the  ex¬ 
treme  complexity  of  the  phenomena.  These  gaps  have  been  tenta¬ 
tively  filled  by  using  approximate  methods  of  calculation,  common 
sense  and  engineering  judgment.  For  example,  a  problem  such  as 
the  penetration  and  propagatio)^  of  a  three-dimensional  bore  in  deep 
estuaries  needs  extensive  theoretical  research  and  development 
which  is  beyond  the  scope  of  Uiis  study  and  report. 

It  is  pointed  out  that  the  general  conclusions  of  this  report 
would  certainly  not  need  to  be  changed  as  a  result  of  further  inves¬ 
tigations,  but  these  conclusions  could  be  refined,  improved  and 
pres ented  with  more  accuracy  subsequent  to  further  theoretical  and 
experimental  investigations. 

It  is  recalled  that  at  the  beginning  ot  the  stiidy  it  was  decided 
that  this  work  would  be  accomplished  in  two  phases:  Phase  I  would 
define  the  problem,  determine  whether  a  problem  actually  existed, 
and  give,  by  use  of  available  techniques  and  formulas,  a  rough 
quantitative  estimate  of  the  amount  of  damage  probable  in  the 
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Chesapeake  Bay  area. 

Phase  11  of  the  study  was  designed  for  improving  the  state 
of  the  art  by  constructing  new  theories.  These  new  theories  have 
been  applied  to  improve  the  accuracy  of  the  final  results  presented 
in  Parts  A  and  B.  Some  of  the  new  theories  have  been  formulated 
only.  Their  application  wovild  require  the  use  of  high  speed  digi¬ 
tal  computers  with  programming  based  on  the  eqviations  which 
have  been  established. 

This  report  is  composed  of  two  main  parts  and  seven  ap¬ 
pendices.  Part  A  is  essentially  nontechnical  and  explains  the 
purpose  of  the  study  and  presents  broad  conclusions  for  the  lower 
Chesapeake  Bay.  Part  B  is  more  technical  and  gives  all  the 
essential  aspects  of  the  problem,  a  brief  presentation  of  the 
methods  of  calculation,  and  their  application  to  the  Norfolk  area. 
The  methods  of  calculation  are  justified  in  indicated  references  or 
in  the  appendices  which  follow.  Some  considerations  on  the 
accuracy  of  these  different  methods  are  also  developed  and  sugges¬ 
tions  given  for  further  investigations.  The  various  appendices  are 
essentially  of  a  theoretical  nature.  Most  of  them  are  original 
contributions  to  the  science  of  hydrodynamics. 

Appendix  I  is  a  synthesis  on  a  critical  literature  survey  in¬ 
cluding  theoretical  and  experimental  information  on  wave  motion 
due  to  a  local  disturbance.  The  variation  of  wave  heights  and 
wave  periods  with  distance  is  also  investigated  quantitatively  for 
Case  III. 
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Appendix  11  reviews  the  now  classical  theory  of  Kranzer 
and  Keller  for  waves  caused  by  explosions  and  its  application  to 
the  problem  under  study.  This  theory  has  been  applied  with  the 
input  data  at  20  miles  from  the  explvosion  as  provided  to  NESCO  by 
United  Reftoarch  Servicoe  for  the  three  cases. 

Appendix  III  presents  the  establishment  of  the  surface  waves 
generated  by  a  disturbance  on  the  sea  bed  in  constant  water  depth. 
This  disturbance  is  defined  as  a  circular  upthrust  varying  ex¬ 
ponentially  with  time  with  an  extremely  rapid  decay  time.  It  is 
an  attempt  at  refinement  of  the  Kranzer  and  Keller  theory  for 
waves  generated  by  underwater  explosions.  Similar  results  for  the 
wave  height  decay  with  distance  arc  presented. 

Appendix  IV  is  also  a  new  approach  for  calculating  the  wave 
motion  generated  by  a  surface  disturbance  of  finite  magnitude  by 
application  of  the  principle  of  superposition  to  the  classical 
Cauchy- Poisson  theory  for  an  infinitely  small  disturbance. 

Appendix  V  deals  with  the  propagation  of  the  wave  on  the  con¬ 
tinental  shelf.  The  shoaling  and  bottom  friction  effects  are  con¬ 
sidered.  Then  a  new  theory  on  ‘'saturated*'  and  "nonsaturated" 
breakers  which  follow  the  breaking  index  ciirve  is  established.  A 
literature  survey  oi\  the  wave  run-up  of  long  waves  on  an  open 
coast  has  been  summarized. 

Appendix  VI  presents  a  new  method  of  computation  of  the  wave 
run-up  based  on  the  method  of  characteristics.  Various  related 
phenomena  have  been  investigated  where  required  for  carrying  out 
successive  step.s  in  the  calculations.  In  particular  the  problem  of 
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the  wave  run-up  on  a  dry  bed  has  been  analy?jed.  This  method 
has  been  applied  to  one  particular  case.  The  systematic  appli¬ 
cation  of  this  method  for  various  bottom  slopes  would  require 
the  use  of  a  high  speed  digital  computer. 

Appendix  VII  presents  a  new  method  of  numerical  compu¬ 
tation  for  the  penetration  of  nonlinear  two-dimensional  waves  into 
bays  and  estuaries.  This  method  would  also  require  use  of  a 
computer  for  its  practical  application  to  the  problem  under  study. 

The  theories  and  methods  developed  in  these  various 
appendices,  and  in  particxilar  in  Appendices  VI  and  VII,  may  have 
many  other  applications  than  this  immediate  problem;  for  example 
tsunami  run-up,  waves  due  to  the  breaking  of  a  dam  and  tidal 
motion,  among  others. 
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where  the  bottom  slope  is  0,01  near  the  shoreline.  This  corres¬ 
ponds  roughly  to  a  twenty-foot  high  breaker.  It  is  also  the  maxi¬ 
mum  possible  average  wave  height  for  the  wave  penetrating  within 
Chesapeake  Bay.  This  maximum  effect  will  occur  from  Cape  Cod 
to  Cape  Hatteras.  However,  the  larger  yield  (Case  III)  located 
further  from  the  coast  will  cause  more  damage  because  its  effects 
will  be  spread  on  a  longer  shoreline,  which  could  even  extend 
from  Canada  to  Florida. 

Considering  the  present  assumed  location  for  the  initial  dis¬ 
turbance,  the  following  effects  can  be  feared. 

Case  I:  The  amount  of  damage  due  to  Case  I  is  very 

limited  and  is  similar  to  the  damage  due  to  a  very  strong  storm. 

The  open  coast  will  be  subjected  to  a  strong  attack  and  the  protecting 
dunes,  which  are  over  15  feet  high  above  the  sea  level,  will  be  sub¬ 
merged.  Inside  Chesapeake  Bay  the  damage  will  be  limited  to 
unprotected  areas  such  as  Old  Point  Comfort  (sec  enclosed  map: 

Fig.  3). 

However,  considering  the  range  of  wave  periods  generated 
by  the  disturbance  and  the  range  of  natural  oscillation  of  the  harbor 
basins,  resonant  motion  will  certainly  occur.  This  motion  will 
have  a  large  amplitude  and  could  cause  considerable  damage  to 
moored  ships.  Hence  this  problem  of  resonance  is  the  most  impor¬ 
tant  factor  in  the  case  of  a  disturbance  of  relatively  small  yield  as 
in  Case  I.  This  problem  can  be  partly  solved  in  some  locations. 

A  breakwater  between  Tanner  Point  and  Craney  Island  can  be 
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considered  as  a  possible  solution  for  the  harbor  basins  in  the 
Elizabeth  River.  An  increase  in  the  length  and  height  of  the 
causeway  between  Old  Point  Comfort  and  Willoughby  Spit  will 
have  a  greater  effect  but  would  probably  be  a  more  expensive 
solution.  However,  this  solution  may  be  worth  consideration  since 
both  .sets  of  harbor  basins  at  Newport  News  and  Norfolk  would  be 
protected.  It  is  evident  that  the  choice  of  such  a  solution  would 
require  further  studies  than  those  presented  in  this  report.  Scale 
model  studies  would  be  highly  desirable. 

Another  .solution  might  be  more  economical.  It  would  con- 
aiat  of  adapting  the  moorings  of  ships  in  such  a  way  that  the  period 
of  oscillation  of  the  mooring  system  does  not  match  the  natural 
period  of  oscillation  of  the  harbor  basin.  This  solution  will  also 
require  a  special  study,  despite  the  fact  that  considerable  work 
has  been  done  on  this  problem  in  the  last  decade. 

Finally,  since  the  waves  would  take  about  one  hour  to  travel 
from  the  entrance  of  Chesapeake  Bay  to  the  harbor  basin  located 
in  the  Norfolk  area,  an  alarm  system  might  be  set  up  in  such  a 
manner  that  certain  precautionary  measures  could  be  taken. 

Cases  II  and  III:  It  has  already  been  pointed  out  that  the 

damage  due  to  these  two  different  yields  located  at  the  same  dis¬ 
tance  from  the  coast  will  be  very  similar  in  the  Norfolk  area. 

Hence  the  results  can  be  presented  togethez’.  This  result  can  be 
considered  as  the  maximum  possible  damage. 
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The  enclosed  map  (Fig.  3)  presents  an  overall  picture  of 
the  areas  subject  to  destruction,  flooding,  and  harbor  resonance. 
First  it  can  be  anticipated  that  the  open  coast  will  be  subjected 
to  such  a  strong  .attack  that  the  protecting  dunes  will  be  submerged 
and  destroyed.  The  shoreline  installations  will  he  destroyed  by  the 
dynamic  effects  of  the  successive  bores  generated  by  the  waves. 

The  flooding  area  will  extend  deeply  within  the  coastal  area  because 
the  elevation  is  generally  low.  Cape  Charles  will  be  half  covered 
by  water;  but  the  most  important  parts  (roads)  which  are  above 
the  EO-foot  contours  will  r  ot  be  submerged. 

The  waves  penetrating  within  Chesapeake  Bay  will  first 
cause  a  slight  rise  in  the  mid  water  level  of  about  one  foot.  But 
this  rise  is  too  small  to  be  a  real  cause  of  concern.  (It  is  even 
smaller  than  the  rise  due  to  storm  surge  which  sometimes  causes 
flooding  in  downtown  Norfolk.)  Hence  the  real  danger  would  result 
from  the  dynamic  effects  of  long  waves.  These  effects  will  be 
particularly  important  from  Cape  Henry  to  Willoughby  Spit.  This 
is  the  first  area  attacked  by  the  waves.  In  this  area  waves  will 
also  pass  over  the  dunes  and  through  Lynnhaven  Inlet  and  Little 
Creek.  The  harbor  installations  near  Little  Creek  will  be  sub¬ 
merged.  The  installations  near  Ocean  View  will  not  only  be  sur- 
merged  but  destroyed  because  of  the  dynamic,  effects  of  the  waves. 

On  the  other  side  of  Hampton  R-oads,  i.  e.  north  and  west  of 
Old  Point  Comfort,  the  dynamic  effects  of  the  waves  arriving 
directly  from  the  entrance  of  Chesapeake  Bay  are  superimposed  on 
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the  waves  reflected  along  the  coast  at'the  south  of  Chesapeake  Bay. 
Hence  a  great  amount  of  destruction  can  be  anticipated  along 
Buckroe  Beach  and  between  Newport  News  and  Hampton. 

Flooding  can  be  anticipated  in  local  places  in  the  downtown 
area  of  Newport  News  and  Norfolk  because  of  resonance  of  harbor 
basins . 

Finally  the  resonance  problem  of  the  harbor  basin  is  similar 
to  the  problem  previously  described  for  Case  I.  It  will  be  greater 
in  amplitude  because  the  incident  wave  height  itself  is  greater. 

The  destructive  effect  of  the  waves  decreases  very  quickly  north 
from  a  line  from  Back  Rover  to  Cape  Charles.  They  can  cause 
some  flooding  in  the  lower  part  of  the  Day. 

The  bridge-tunnel  between  Chesapeake  Beach  and  Fisherman's 
Island  will  not  change  the  overall  results  presented  in  this  report 
but  this  bridge- tunnel  can  be  undermined  by  erosion  and  subjected 
to  destruction  because  of  wave  forces. 

Now  the  military  installations  in  the  Norfolk  area  will  be 
considered  more  specifically.  Most  of  the  military  installations 
are  subjected  to  wave  effect  because  they  are  located  in  the  lowest 
areas.  The  probably  maximum  damage  can  be  summed  up  as  follows 

1 .  Langley  Air  Force  Base 

Local  flooding  around  shoreline.  No  dynamic  destruc¬ 
tion. 
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Z.  Fort  Monroe 

Dynamic  destruction  along  shoreline  accompanied 
by  extensive  flooding  throughout. 

U>  S.  Military  Reservation  at  Newport  News 
Some  flooding  near  the  waterline. 

4.  U.  S.  Naval  Reservation  at  Willoughby 

Some  dynamic  damage  at  SewelFs  Point.  Extensive 
flooding.  Seiche  motion  of  high  amplitude  in  harbor 
basins . 

5 .  U.  S.  Naval  Reservation  at  Craney  Island 
Extensive  flooding. 

6.  Naval  Hospital  at  Hospital  Point 

Little  or  no  damage, 

7.  Norfolk  Naval  Shipyard 

Resonance  of  harbor  basins  causing  seiche  motion 
and  possible  local  flooding  due  to  the  resonance. 

8.  U,  S.  Military  Reservation  near  Paradise  Creek 
Resonance  of  harbor  bavsins  causing  seiche  motion 
and  poi,sible  local  flooding  due  to  the  resonance. 

9 .  U.  S.  Naval  Re.servation  at  Little  Creek 
Dynamic  destruction  along  the  shoreline.  Extensive 
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2.  Fort  Monroe 

Dynamic  destruction  along  shoreline  accompanied 
by  extensive  flooding  throughout. 

3.  U.  S.  Military  Reservation  at  Newport  News 
Some  flooding  near  the  waterline. 

4.  U,  S.  Naval  Reservation  at  Willoughby 

Some  dynamic  damage  at  Sewell's  Point.  Extensive 
flooding.  Seiche  motion  of  high  amplitude  in  harbor 
basins . 

5 .  U.  S.  Naval  Reservation  at  Craney  Island 
Extensive  flooding. 

6.  Naval  Hospital  at  Hospital  Point 

Little  or  no  damage. 

7 .  Norfolk  Naval  Shipyard 

Resonance  of  harbor  basins  causing  seiche  motion 
and  possible  local  flooding  due  to  the  resonance. 

8.  U.  S.  Military  Reservation  near  Paradise  Creek 
Resonance  of  harbor  basins  causing  seiche  motion 
and  possible  local  flooding  due  to  the  resnnanr-R 

9.  U.  S.  Naval  Reservation  at  Little  Creek 
Dynamic  destruction  along  the  shoreline.  Extensive 
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flooding  around  Little  Craek,  resonance  of  harbor, 
basin,  the  navigation  channel  filled  by  sand. 

10.  U.  S,  Military  Reservation  at  Camp  Pendleton 
Little  or  no  local  flooding.  No  damage. 

11.  U.  S.  Naval  Reservation  South  of  Virginia  Beach 
Complete  erosion  of  the  dunes.  Complete  dynamic 
destruction  along  the  shoreline.  Complete  flooding. 

12.  U.  S.  Military  Reservation  at  Fort  Story 
Dynamic  destruction  along  the  shoreline.  Partial 
flooding. 

1^*  Cape  Charles  Air  Force  Base 

Some  dynamic  destruction  and  extensive  flooding. 

Any  possible  solution  to  the  problem  on  the  open  coast 
would  be  unrealistic.  Engineering  solutions  could  only  be  con¬ 
sidered  for  local  protection.  Some  solutions  have  already  been 
outlined  in  Case  I.  The  solutions  for  protection  against  conditions 
arising  from  Cases  II  and  III  will  be  similar  and  also  would  require 
further  investigations,  involving  scale  model  studios. 

Finally*  it  may  be  pointed  out  that  the  results  presented  here 
ha.ve  been  obtained  first  by  application  of  existing  wave  theories. 

The  difficulty  of  the  problem  is  such  that  a  number  of  new  theories 
were  necessary  in  order  to  improve  the  state  of  the  art.  Despite 
these  new  theories,  many  points  rely  on  engineering  judgment  rather 
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than  scientific  calculation.  The  results  obtained  have  indicated 
the  extent  of  the  problem.  They  can  be  subjected  to  slight  change 
after  further  analysis.  There  is  no  doubt  that  more  theoretical 
and  experimental  studies  in  the  field  or  in  a  laboratory  would  give 
considerable  complementary  information  of  great  value. 

The  conclusions  can  be  summarized  in  a  few  words.  The 
destructive  action  of  long  waves  generated  by  a  nuclear  explosion 
off  the  Atlantic  coast  would  be  lees  serious  than  may  be  expected, 
a  priori,  because  of  the  natural  protection  afforded  by  the  con¬ 
tinental  shelf.  On  the  other  hand,  the  possible  destruction  would 
be  much  more  extensive  than  that  due  to  the  strongest  observed 
storm.  Hence  the  destruction  may  be  sufficient  to  be  a  cause  of 
concern. 

Finally,  It  is  pointed  out  that  the  protective  action  of  the 
continental  shelf  arises  as  a  consequence  of  its  length  and  slope. 
Hence  the  conclusions  presented  in  this  report  cannot  be  extrapo¬ 
lated  to  the  Pacific  coast  because  there  the  continental  shelf  is 
much  shorter  and  steeper. 
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PART  B 


DESCRIPTION  OF  PHYSICAL  PHENOMENA,  METHODS 
OF  COMPUTATION,  RESULTS  AND  RECOMMENDATIONS 
FOR  FURTHER  INVESTIGATIONS 


By 

B.  Le  Mehaute 
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CHAPTER  in 

A  DESCRIPTION  OF  THE  PHYSICAL  PHENOMENA 

The  detailed  description  of  the  wave  motion  at  the  location  of 
the  disturbance  is  not  within  the  scope  of  this  study.  However,  for 
a  better  understanding  of  the  subsequent  phenomena  which  concern 
this  report,  a  brief  description  is  presented. 

Consider  the  case  of  a  sudden  disturbance,  such  as  that  due 
to  the  explosion  of  a  nuclear  device  below  the  sea  surface  in  deep 
water.  From  a  hydrodynamic  point  of  view  this  disturbance  could 
be  considered  as  due  to  successive  expansions  and  contractions  of 
a  huge  cavity.  A  hot  gas  and  steam  bubble  is  quickly  created  which 
initiates  shock  waves.  Upon  arrival  at  the  free  surface  these  shock 
waves  cause  the  formation  of  a  spray  dome  which,  when  falling, 
causes  a  '^base  surge.  “  Then  the  hot  gas  and  steam  bubble  burets 
through  the  dome  causing  considerable  water  movement. 

Both  theory  and  experiment  yield  methods  for  the  determina¬ 
tion  of  the  wave  motion  far  from  the  location  of  the  disturbance. 

All  the  theories  are  based  on  a  relatively  simple  mathematical 
model  assuming  the  depth  to  be  either  constant  or  infinite.  This 
assumption  of  infinite  depth  is  reasonably  justified  for  the  wave 
somewhere  near  the  disturbance. 

In  spite  of  the  complexity  of  the  real  motion  and  great  sim¬ 
plification  of  the  mathematical  models,  the  essential  features  of 
the  waves  given  by  theory  are  relatively  well  verified  by  teats, 
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These  waves  are  characterized  by  relatively  long  periods,  usually 
ranging  from  30  seconds  to  a  few  minutes.  They  are  slightly 
shorter  but  much  steeper  than  tsunami  waves  due  to  seismic  dis¬ 
turbances.  However,  like  tsunami  waves,  the  length  between 
crests  ie  so  great  that  their  steepness  in  deep  water  is  often  in¬ 
significant  and  they  are  seldom  visible  some  distance  from  their 
source.  The  waves  are  circular  and  travel  radially  from  the 
origin  in  a  manner  similar  to  waves  caused  by  a  stone  dropped  in 
a  pond.  This  results  in  the  wave  energy  being  spread  over  a  circle 
of  increasing  radius,  and  the  wave  height  decreases  rapidly. 

Such  waves  are  not  periodic  in  the  mathematical  sense: 
their  "pseudo  period"  varies  with  time  and  distance.  Due  to  the 
increase  in  wave  period  with  distance,  the  "wave  length"  (or  more 
exactly  the  distance  between  wave  crests)  and  the  group  velocity 
{or  velocity  of  propagation  of  energy)  also  increases  with  distance. 

In  other  words,  the  "wave  lengths"  are  "stretched"  radially.  This 
stretching  effect  is  particularly  important  in  deep  water.  Stretch¬ 
ing  and  spreading  effects  cause  the  average  wave  height  to  decrease 
very  quickly  with  the  distance  from  the  origin  of  the  disturbance. 

The  Btrerching  effect  becomes  less  important  in  shallow  water  and 
becomes  negligible  in  very  shallow  water,  such  as  on  the  continen¬ 
tal  shelf.  In  very  shallow  water  the  wave  decay  is  due  only  to  the 
"spreading"  effect  and  the  wave  height  decreases  in  a  manner 
similar  to  an  ordinary  periodic  gravity  wave  in  constant  water  depth. 

The  wave  height  and  wave  period  also  vary  from  one  wave  to 
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another.  The  wave  motion  appears  as  a  succession  of  wave  train 
presenting  a  "beating  effect."  The  number  of  waves  in  a  wave 
train  also  increases  with  time  and  distance.  The  number  of  maxi¬ 
mum  waves  also  increases  with  distance  from  the  origin.  In 
practice,  the  wave  wliicii  is  the  liigUesl  in  the  wave  train  at  any 
given  instant  retau^s  the  same  period. 

The  first  wave  train  ia  the  most  important  because  it  con¬ 
tains  the  most  energy  and  the  highest  waves.  The  general  ten¬ 
dency  at  a  given  point  ia  that  the  wave  period,  height  of  the  envelope 
of  wave  trains,  number  of  waves  within  a  wave  train,  and  length 
of  the  wave  train  decrease  with  time.  Also,  the  wave  height  at  a 
given  point  and  a  given  time  depends  upon  the  yield  of  the  distur¬ 
bance,  while  the  wave  period  depends  upon  time  and  distance  only. 

Because  the  sea  bottom  is  not  horiz.ontal>  these  wave  trains 
are  theoretically  influenced  by  a  "shoaling  effect."  The  shoaling 
effect,  due  to  changes  in  depth,  is  added  to  this  spreading  and 
stretching  effect  on  a  bottom  of  varying  depth.  The  decrease  in 
depth  causes  an  increase  in  wave  height  which  partly  compensates 
for  the  wave  decay  due  to  spreading  and  stretching  and  becomes 
particularly  important  upon  arrival  on  the  continental  slope.  It 
is  interesting  to  note  that  a  shoaling  effect  also  exists,  in  a  sense, 
for  a  constant  depth.  Since  the  "wave  length"  increases  with 
distance  from  the  origin,  a  wave  motiuit  near  tiiC  source  in  deep 
water  travels  on  a  continually  smaller  "relative  depth"  and  will 
even  tend  to  travel  in  relatively  shallow  water  at  an  infinite 
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slope.  However,  in  practice  the  valife.s  found  from  both  as.sump- 
tions  are  smaller  than  double  the  deep  water  wave  height  for  this 
particular  problem.  The  truth  lies  between  the  two  extremes  of 
reflection  on  a  vertical  wall  and  total  energy  transmission.  No 
reliable  theories  are  available  to  give  the  exact  feature  of  the 
change  in  wave  height.  However,  an  attempt  is  made  to  give  a 
reasonable  answer  to  this  problem.  It  has  been  found  that  the 
two  extremes  corresponding  to  either  of  the  two  assumptions  above 
give  results  close  enough  to  permit  a  very  good  estimation  of  the 
effect  of  the  continental  slope. 

Over  the  continental  shelf  the  waves  can  first  be  considered 
as  periodic  waves.  The  wave  height  tends  to  increase  due  to 
shoaling.  At  the  same  time,  the  bottom  friction  becomes  more 
and  more  important  and  causes  damping,  partially  compensating 
for  the  increase  of  wave  height  due  to  shoaling.  For  the  highest 
waves  the  bottom  friction  ie  so  important  that  the  wave  height  can 
even  decrease  despite  the  shoaling  effect.  At  the  same  time  that 
the  wave  height  tends  to  increase,  the  "wave  length"  or  distance 
between  two  successive  wave  crests  is  decreasing.  As  a  conse¬ 
quence,  the  wave  becomes  steeper  and  steeper. 

These  so-called  pseudo- periodic  waves,  which  have  an  almost 
sinusoidal  profile  in  deep  water,  then  appear  as  a  succession  of 
solitary  waves  which  are  characterized  by  high  peaks  of  water  be¬ 
tween  very  long,  flat  troughs.  These  solitary  waves  are  not  oscil¬ 
latory  waves,  for  which  the  flow  is  periodic  and  the  total  discharge 
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over  the  period  of  tlie  wave  13  theoretically  zero.  They  are 
acUially  translatory  waves,  characterized  by  a  transportation  of 
water  in  the  wave  direction  with  no  return  flow.  Actually  the 
transport  of  water  with  the  crest  of  a  solitary  wave  causes  a  rise 
in  the  mean  water  level  near  the  shoreline  resulting  in  a  slow 
return  flow  in  the  succeeding  long,  flat  trough.  Despite  the  re¬ 
turn  flow  the  wave  profile  is  actually  very  close  to  that  of  a 
solitary  wave  when  the  ratio  of  the  depth  and  the  distance  between 
crests  becomes  small  and  the  bottom  slope  is  very  gentle,  condi¬ 
tions  encountered  on  the  continental  shelf  of  the  Atlantic  coast. 

The  wave  height  in  a  solitary  wave  also  increases  when 
the  depth  decreases  but  theoretically  this  change  is  much  greater 
than  for  a  periodic  wave.  The  increase  in  wave  height  depends 
upon  the  slope  and  the  relative  depth.  However,  a  reliable 
mathematical  law  to  predict  the  increase  accurately  is  unknown. 
The  theory  based  on  the  conservation  of  energy  in  a  solitary  wave 
is  not  verified  by  experiment.  Long  waves  are  partially  reflected 
on  a  very  gentle  slope.  At  present  the  best  method  available,  as 
a  preliminary  approximation,  is  to  extrapolate  available  experi¬ 
mental  data.  However,  it  must  be  borne  in  mind  that  such  extra¬ 
polation  is  always  subject  to  doubt  or  at  least  to  qualitative 
interpretation. 

As  the  solitary  wave  becomes  steeper,  its  celerity  becomes 
strongly  influenced  by  its  amplitude.  As  the  depth  decreases  the 
particle  velocity  at  the  crest  tends  to  become  greater  than  the  wave 
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energy  is  dissipated  by  a  spilling  breaker  before  reaching  the 
shoreline. 

The  greater  the  wave  height  in  deep  water,  the  sooner  the 
wave  will  break.  Hence,  the  amount  of  damage  caused  by  this 
wave  is  almost  independent  of  the  initial  height  of  the  wave  in 
deep  water  once  this  wave  height  exceeds  a  critical  value.  This 
critical  value  of  wave  height  is  characterized  by  a  slope  steep 
enough  for  the  breaker  to  become  saturated  or  fully  developed 
bore.  Such  a  slope  exists  only  near  the  shoreline  in  relatively 
very  shallow  water.  {It  can  also  exist  on  the  continental  elope  but 
either  the  water  is  too  deep  to  cause  any  type  of  breaking  or,  even 
if  breaking  occurs  for  huge  waves  over  the  continental  slope,  the 
saturated  breaker  again  becomes  a  nonssturated  breaker  upon 
reaching  the  gentle  slope  of  the  continental  shelf.  Nonsaturated 
breakers  and  even  nonbreaking  waves  can  again  reform  from 
this  breaker.  )  Hence  the  continental  shelf  off  the  Atlantic  coast 
acts  as  a  natural  protection.  This  natural  protection  would  be 
much  less  efficient  on  the  Pacific  coast  because  of  the  steeper 
elope  of  the  continental  shelf. 

The  maximum  damage  can  occur  either  from  a  breaker  initially 
spilling  all  along  the  continental  shelf  or  from  a  wave  of  initially 
smaller  amplitude  in  deep  water  which  reaches  its  breaking  condi¬ 
tion  at  a.  critical  slope  where  the  breakers  are  saturated.  There¬ 
fore,  despite  the  natural  protection  of  the  continental  shelf,  the  slope 
near  the  shoreline  becomes  steep  and  the  breaker  is  satura-ted  at  a 
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depth  where  the  wave  is  high  enough  (taay  2,0  feet)  to  be  a  cause 
of  great  concern,  in  which  case  the  wave  is  transformed  into  a 
translatory  wave  such  as  a  tidal  bore. 

No  simple  relationship  exists  to  relate  the  height  of  a 
breaking  solitary  wave  to  Uie  height  at  the  front  of  the  tidal  bore 
issued  from  this  wave.  However,  they  are  of  the  same  order  of 
magnitude  near  the  point  where  the  breaker  becomes  saturated. 

The  tidal  bore  moves  at  a  speed  which  depends  upon  the  height 
of  the  wave  front  and  the  depth. 

As  in  any  shock  wave,  a  great  amount  of  energy  will  be 
dissipated  near  the  front  of  the  bore  by  sudden  changes  in  depth. 

The  bottom  friction  becomes  relatively  less  important  than  the 
internal  friction  within  the  bore.  The  front  of  the  tidal  bore 
could  be  considered  as  a  moving  hydraulic  jump  with  the  velocity 
of  the  upstream  flow  equal  to  zero.  The  relative  amount  of  wave 
energy  dissipated  increases  as  the  depth  decreases.  The  bottom 
friction  becomes  more  and  more  important  as  the  wave  nears 
the  shoreline. 

From  a  hydrodynamic  point  of  view  the  bore  considered  as 
a  shock  wave  does  not  exist  any  more  when  climbing  on  a  dry  bed. 

It  is  a  "rarefaction  wave"  or  edge  of  water  which  moves  upwards. 
But  the  bottom  friction  forces  tend  to  be  infinite  at  the  extreme  edge 
of  the  wave.  Hence,  in  practice  this  edge  is  cut  short  and  has  the 
physical  appearance  of  a  bore. 

However,  the  rise  of  the  wave  on  a  dry  bed  cannot  be  treated 
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from  incident  waves.)  According  to  the  linear  theory,  along  the 
separation  line  the  wave  amplitude  has  half  of  the  value  at  the  en¬ 
trance.  In  the  shaded  areas  the  breaker  can  generate  a  non¬ 
breaking  wave  of  much  smaller  wave  height.  These  new  waves 
will  break  when  they  arrive  on  the  ahoreline  generating  a  bore. 

This  bore,  while  being  much  smaller  than  on  the  open  coast  and 
in  unprotected  areas,  wovild  still  be  large  enough  in  many  places  to 
cause  considerable  flooding  and  damage  on  a  flat  shoreline. 

Finally,  the  harbor  basins  could  be  subjected  to  resonance 
effects.  Available  theories  permit  the  calculation  of  the  periods  of 
resonance,  and  even  amplification  factors,  with  relatively  good 
accuracy.  Most  harbor  basins  are  similar  to  an  open  acoustic 
pipe.  They  can  be  subjected  to  a  resonance  phenorr.enon  when  the 
period  of  the  waves  at  the  entrance  matches  one  of  the  natural 
periods  of  free  oscillation  of  the  basin.  This  occurs  when  the 
length  of  the  open  basin  is  close  to  a  quarter  of  the  wave  length 
(or  3/4,  5/4,  .  .  .  for  harmonics).  Considering  the  usual  siae  and 
depth  of  the  harbor  basins,  it  has  been  found  tl^at  most  of  them  will 
be  strongly  excited  by  the  range  of  periods  of  the  expected  incident 
waves . 

Since  the  incident  waves  are  not  periodic  in  the  mathematical 
sense,  the  problem  of  resonance  in  a  harbor  must  be  considered 
from  a  spectrum  point  of  view.  The  incident  waves  can  be  considered 
as  defined  by  a  wave  spectrum.  The  time  of  response  of  an  open 
basin  to  this  incident  spectrum  is  very  short  (three  or  four  wave 
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periods)  and  at  resonance  the  amplitiide  is  great  (at  least  four 
times  the  incident  wave  height).  This  can  be  explained  from  the 
fact  that  a  great  amount  of  wave  energy  penetrates  through  the 
entrance,  establishing  the  resonance  condition  very  quickly. 

Hence,  even  in  the  areas  protected  from  destruction  by  the 
dynamics  of  the  waves  and  from  gradual  flooding,  this  problem 
must  also  be  a  cause  of  concern.  This  is  due  to  the  fact  that 
moored  ships  are  then  subjected  to  a  seiche  motion  of  large  am¬ 
plitude.  They  may  act  in  a  similar  manner  to  a  body  between  two 
springs.  In  these  cases  also  the  period  of  resonance  of  the  basin 
is  close  to  the  period  of  resonance  of  large  boats  moored  along 
the  quays.  Consequently  the  moorings  can  be  broken  and  the 
ships  thrown  bodily  against  the  quays  or  each  other. 

Now  that  the  main  physical  phenomena  have  been  qualita¬ 
tively  described,  the  methods  of  calculation  will’ be  briefly  sum¬ 
marized  in  the  following  chapter.  Their  justification  will  be  found 
in  various  appendices  or  listed  references. 
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convective  inertia  term,  which  is  usulally  neglected  in  theories, 
becomes  important  for  steep  waves.  Hence,  it  can  be  expected 
that  the  wave  decay  near  the  origin  is  larger  than  the  decay  given 
by  theories,  namely  H  CC  R  ^  for  the  maximum  of  the  wave  en¬ 
velope  {H  is  the  wave  height  and  R  the  distance  from  the  origin 
of  the  disturbance),  {See  Appendices  i,  II  and  III.) 

The  shoaling  effect  due  to  the  change  of  depth  is  added  to 
the  shoaling  effect  due  to  the  increase  of  wave  length. 

It  may  be  seen  then  that  as  the  waves  spread  out  and  their 
periods  increase  with  distance,  the  leading  waves  which  are  the 
longest  begin  to  feel  the  bottom  first  and  could  be  subject  to  shoal¬ 
ing  effects  before  even  reaching  the  continental  slope.  This  shoal¬ 
ing  effect  due  to  a  change  of  relative  depth  can  be  estimated  at 
least  in  a  first  r  nproxiniation. 

Before  studying  the  shoaling  effect  for  the  first  waves  before 
the  continental  slope,  some  considerations  must  be  given  to  the 
variation  of  the  ’’wave  period”  with  distance.  It  has  been  seen  that 
this  kind  of  motion  is  nonperiodicai  in  the  mathematical  sense. 
However,  a  ’’pseudo”  wave  period  can  be  defined.  This  wave  period 
is  defined  as  the  period  of  a  periodic  wave  having  the  same  distance 
between  wave  crests.  Then  it  is  found  that  this  wave  period  varies 
with  distance.  The  law  of  variation  is  actually  very  complex  (see 
Appendices  I  and  II).  It  varies  both  Iroxn  the  first  wave  to  the  last 
wave  and  with  distance.  As  a  first  approximation  the  wave  period 
for  a  given  wave  may  be  considered  as  increasing  with  the  cube  root 
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It  is  recalled  that  the  group  velocity  V  is  given  by  the 
formula 


V 


C, 

2 


4  ^  d 

L 

sinh 


where  L  ie  the  wave  length  (L  s:  - 

A.  tr 

the  depth  and  C  the  wave  celerity  (G  =  Hence 

at  = 


tanh  m| 

2  m,  d, 

If  ' 

sinh  2m,  d. 

tonh  ntg  dg 

r  ,  2m,  - 

sinh  2mg  dg 

. 

{m  = 

The  cumulative  correcting  factor  to  be  applied  to  the  results 
obtained  from  the  Kranzer  and  Keller  theory  (Appendix  11)  Is  given 

n 

by  their  product;  J"’  =  j  j  .  A  P 

It  can  be  seen  that  due  to  the  change  of  depth,  the  wave  height 
of  the  first  waves  must  be  increased  by  a  small  coefficient.  This 
correction  is  negligible  for  the  highest  waves  which  are  always  in 
deep  water  (in  the  mathematical  sense)  before  reaching  the  continen¬ 
tal  slope.  The  same  method  may  be  applied  to  the  continental  slope 
and  will  give  a  much  higher  coefficient  for  all  waves.  This  is  the 
subject  of  the  following  paragraph. 


Z ,  The  Evolution  of  the  Waves  on  the  Ccntinental  Slope 

All  the  wave  energy  is  not  transmitted  over  the  continental 
shelf.  A  part  of  it  is  reflected  towards  the  sea  by  the  continental 
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slope.  This  problem  is  considered  ip  this  section. 

First,  the  variation  of  the  wave  period  over  the  continental 
slope  is  negligible  and  T  can  be  considered  as  a  constant  for  a 
given  wave.  The  continental  slope  on  the  Atlantic  coast  extends 
from  a  mean  depth  of  100  fathoms  to  a  depth  of  1000  fathoms 
over  an  average  horizontal  distance  of  10  miles,  giving  a  slope 
of  i/100.  This  slope  may  be  steeper  in  local  areas  but  it  ie 
never  steeper  than  1/10. 

It  has  been  demonstrated  that  when  L  »  2  ir  a  ,  where 

a  is  the  luorizontai  distance  (presently  about  10  miles),  the  con¬ 
tinental  slope  may  be  considered  as  a  vertical  wall.  (See  Yoshida 
(1950)  )  This  condition  is  not  satisfied  lor  this  case  but  con¬ 
sidering  the  continental  slope  as  a  vertical  wall  presents  the 
advantage  of  having  a  limit.  The  true  fact  lies  between  this  limit 
obtained  by  considering  that  all  energy  is  transmitted  over  the 
continental  shelf  following  the  method  previously  presented.  It 
will  be  seen  that  (I)  these  two  limits  yield  relatively  close  solu¬ 
tions,  and  (2)  any  error  in  the  deep  water  wave  problem  does  not 
change  the  final  conclusions  because  of  the  regulating  effect  of  the 
continental  shelf. 

Considering  the  continental  slope  as  a  vertical  wall,  the  co¬ 
efficient  of  transmission  a  is:  (See  Le  Mehaute,  (1961)  ) 


a 


+  A-r 


U) 
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a  is  the  ratio  of  the  wave  height  ovpr  the  continental  shelf  and 
the  wave  height  in  deep  water. 


A 


1  + 


I  + 


dg 

sintranT^T^ 


sinh2m|  d. 


(2) 


(Index  1  refers  to  deep  water;  index  2  refers  to  the  continental 
shelf. )  It  may  be  noted  that  in  the  case  of  reUtively  shallow 
water,  A  =  1  and 


a 


(3) 


a  formula  also  proposed  by  Poincare  (Ed,  1931),  Lamb  (1932), 
Cochrane  and  Arthur  (1948)  and  others  for  tidal  waves  or  very 
long  waves . 

Considering  the  usual  approximations,  the  coefficient  of 
transmission  for  the  first  waves  of  the  wave  train  is: 


=  2 


i  + 


[gdj]  '/a 


'i  ' 
2r- 


2m.  d. 


sinh  2m, 


(4) 
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For  the  highest  wave  which  has  a  shorter  period,  this  simplifies 
to: 


f 

Zm,  d,  1 

1  +  —  ■—  f 

■ 

I  ^  \ 

'  tanh  rrig  dg 

1 

sinh  2mj  dg  J 

In  any  case,  the  general  formula  (1)  always  applies. 


The  spreading  effect  will  be  taken  into  account  by  multiply- 
ing  the  above  values  for  a  W  \p — )  for  a  periodic  wave, 


and  being  the  distance  from  the  origin  at  the  toe  of  the 

continental  slope  and  on  the  edge  of  the  continental  shelf,  respec~ 
tively.  The  problem  remains  to  choose  values  of  and 
because  the  variation  of  bottom  slope  is  progressive.  This  point 
will  be  fully  discussed  in  the  following  chapter. 

The  continental  slope  is  connected  to  the  continental  shelf 
by  a  gentle  bend  from  a  depth  of  100  fathoms  to  a  depth  of  20 
fathoms.  Over  this  part  the  Green  Law  applies  {see  Appendix  V): 


(  "i 


TV  ^ 

which  can  be  refined  by  application  of 


1/4 


dj  +  Hi 
^12  + 


1/4 


which  itself  can  also  be  corrected  for  friction  lorces. 

Passing  over  the  continental  slope,  the  waves  arrive  on  the 
continental  shelf. 
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3 .  The  Evolution  of  the  Waves  on  the  Continental  Shelf  Before 

Breaking  (Appendix  V) 

The  conHnental  shelf  of  the  Atlantic  coast  of  the  United  States 
is  characterized  by  a  slope  so  gentle  that  it  can  almost  be  con¬ 
sidered  horizontal.  The  slope  above  the  fathom  contour  varies 
from  2/10,000  to  5/10,000,  becoming  steeper  near  the  shoreline, 
and  ending  with  a  beach  slope  between  1/30  and  1/10.  The  con¬ 
tinental  shelf  extends  over  a  length  ranging  from  50  to  100  miles 
or  more. 

The  periodic  wave  theories  and  solitary  wave  theory  can  be 
used  as  sound  guides  for  studying  wave  motion  over  the  continental 
shelf.  At  first  the  most  convenient  theory  for  treating  the  problem 
at  hand  must  be  chosen.  The  chart  shown  in  Fig,  1-4  of  Appendix  I 
indicates  the  range  of  validity  of  the  various  theories  used  in 
practice.  However,  according  to  some  experimental  results 
presented  in  Fig.  V-2  of  Appendix  V,  it  seems  that  the  Airy  theory 
extends  up  to  a  point  where  -j-  =1.4  where  d^  is  the  breaking 

depth.  It  is  this  last  criterion  which  will  be  adopted.  (See  Munk 
(1949)  i 

Now  the  Airy  law  has  been  modified  to  take  nonlinear  con¬ 
vective  inertia  and  bottom  friction  into  account.  Then  it  is  demon¬ 
strated  (see  Appendix  V)  that  the  forrnvilas  to  be  used  are: 


H,  =  H, 


/d,  +  H,y4  j 

1  _ 

fH,  AR  ' 

/ 

3ir 

1 

i/z 
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where  A  R  is  the  horizontal  dietancp  between  two  cross  sections 
defined  by  the  distance  and  from  the  origin  (AR  “  -  R^) 

depth  and  d.^,  and  f  ie  the  friction  factor  dimensionless. 

If  the  wave  period  is  long  enough,  the  wave  motion  can  be 
considered  as  a  succession  of  steady  motions  as  far  as  the  vertical 
velocity  distribution,  shearing  stress  and  friction  forces  are  con" 
cerned.  In  that  case,  f  =  — where  la  the  Cheay  coeffi- 

^1/6''*’'' 

cient  and  Cj^  ®  1.486  -—-j —  where  n  is  the  Manning  coefficient( 
and  d  is  measured  in  feet.  With  n  =  0.  02,  this  yields; 


2 

3tr 


A  H  =  5  miles  S  30,  000  feet  is  a  convenient  interval  for 
calculating  the  wave  height  step  by  step  over  the  continental 
shelf.  For  shallower  water  A  R  must  be  chosen  as  a  smaller 
value  since  the  formula  for  the  damping  of  the  wave  assumes  that 

^  ^  «  TMinr 

Now,  if  v  ‘-  <  1.4,  the  following  formula  based  on  the 

solitary  wave  theory  can  be  applied:  (See  Appendix  V  and  Ippen 
and  Kulin  (1954)  ) 


up  to  the  point  where  H  reaches  a  value  ^  0.  78  d.  Then  the 

wave  breaks . 
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Such  is  rarely  the  case.  The  'yave  arrives  at  the  entrance 
to  the  bay  and  estuary  with  a  wave  height  equal  to  0.78  d.  This 
result  holds  true  whatever  the  deep  water  wave  height,  provided 
it  is  high  enough  to  reach  the  breaking  index  curve  at  the  con¬ 
sidered  depth. 

The  maximum  run-up  can  be  calculated  by  considering  the 
breaking  wave  height  (H^  n  0, 78  d)  at  the  location  where  the  breaker 
becomes  saturated,  i.  e.  where  the  slope  becomes  steeper  than 
0.  01.  This  holds  true  provided  the  slope  maintains  this  constant 
value  or  becomes  even  steeper.  It  cannot  be  true  if  the  slope  be¬ 
comes  less  than  0.01.  The  breakers  reform  into  nonsaturated 
breakers  over  this  slope.  A  saturated  breaker  on  the  continental 
slope,  if  any,  wilt  become  nonsatur.ated  over  the  continental 
shelf  and  finally  reach  the  breaking  index  curve  H  =  0.78  d  as 
a  smaller  wave. 

It  has  been  found  that  the  run-up  R  of  a  long  wave  over  a  I/IO 
siope  ib  given  by  s  2.  82,  H  being  the  wave  height  where  the 
bottom  slope  is  0,01,  i.c,  0.78  dj^  [see  Appendices  V  and  VI). 
Other  formulas  based  on  experimental  observations  for  slopes  of 
1/30  and  1/60  are  also  reviewed  in  Appendix  V  and  indicate  similar 
results . 

5.  Penetration  of  the  Waves  Into  Estuaries 
(Appendix  VII) 

The  penetration  of  waves  into  the  lower  Chesapeake  Bay  may 
be  considered  as  the  penetration  of  a  wave  Into  a  basin,  as  shown 
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by  Fig.  4.  The  wave  diffracts  north'xyard  and  quickly  becomes 
negligible.  According  to  the  linear  theory  of  wave  diffraction, 
the  wave  height  along  AB  is  i/2  of  the  wave  height  at  the  entrance 
and  decreases  very  quickly  northwards. 


FIGURE  4 

SCHEME  OF  LOWER  CHESAPEAKE  BAY 

Although  the  wave  motion  under  consideration  is  highly  non¬ 
linear,  such  results  remain  qualitati veiy  true,  but  nonlinear  long 
waves  and  bores  will  be  expected  to  give  more  wave  agitation  than 
is  indicated  by  the  linear  theory. 

As  a  first  approximation,  the  side  CD  can  be  considered  as 
a  vertical  wall.  But  the  usual  variations  of  wave  height  given  by 
the  theory  on  the  side  BD,  which  theoretically  range  from  1,175 
to  0.  88,  cannot  be  expected  to  be  observed  in  this  case  because  of 
the  complexity  of  the  sea  bottom. 

The  wave  calculation  from  AC  to  BD  should  instead  take 
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Then  the  maximum  wave  height  (and  limit  steepness)  of  short 
crested  waves  is  much  higher  than  in  the  case  of  long  crested 
waves.  A  formula  has  been  proposed  for  the  limit  steepness  of 
auperim posed  waves,  whatever  their  relative  direction.  It  applies 
for  irregular  waves,  ciapotis,  or  short  crested  waves.  (See 
Le  Mehaute  (196i)  ).  For  N  waves  this  formula  is.: 


Z 


tanh 


gird 

L 


In  the  case  of  long  waves  this  formula  gives: 

H  =  2  H.  =  1.256  d 
max  b 

In  the  second  case,  where  the  energy  diverges,  the  wave 
celerity  is  (d  +  H)  where  H  can  be  calculated  by  applying 

die  principle  of  conservation  of  energy  between  orthogonals  from 
the  considered  location  and  the  previous  location  where  H  departs 
from  the  breaking  index  condition,  i.  e. 

In  this  formula  b^,  and  are  the  distances  between  orthogonals. 

In  fact  the  principle  of  conservation  of  energy  between  ortho- 
gonals  must  also  be  corrected  empirically  to  take  into  account 
diffraction  effects.  Hence  in  "shaded'*  areas  the  wave  breaks  only 
near  the  shoreline. 

Now  consider  more  carefully  the  side  DC,  which  In  fact  is 
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curved.  A  part  of  the  wave  energy  will  be  dissipated  by  lateral 
breaking  and  penetration  over  land  {see  Chen  (1945)  ).  A  part 
will  be  reflected  towards  side  BD.  This  reflected  part  will  be 
superimposed  on  the  wave  arriving  directly  from  AC,  But  also, 
along  the  side  DC  the  so-called  phenomenon  of  a  stem  or  edge 
wave  may  be  expected  because  the  angle  of  incidence  is  small. 

This  phenomenon  of  stem  wave  partly  reduces  the  loss  of  energy 
from  C  to  D. 

The  problem  of  wave  run-up  along  the  side  BD  is  theo¬ 
retically  very  similar  to  the  problem  of  the  wave  run-up  on  the 
open  coast.  In  practice  the  problem  is  quantitatively  different 
because  the  depth  near  the  shoreline  is  emailer.  Inside  the 
smaller  estuaries,  near  Norfolk,  the  problem  is  that  of  succes¬ 
sive  diffraction  with  multiple  reflection  and  damping  due  to 
bottom  friction. 

The  ratio  of  width  of  the  channel  to  the  wave  length  is  always 
too  large  to  consider  the  motion  as  similar  to  that  of  a  tidal  wave. 
Hence  the  diffraction  phenomenon  assumes  great  importance.  Un¬ 
fortunately  very  little  is  known  about  the  diffraction  of  long  waves. 
The  method  proposed  in  Appendix  VII  can  be  applied  to  solve  this 
problem . 

Another  problem  which  can  be  a  cause  of  concern  is  the 
elevation  of  the  jtij*::!  •,  t  .cr  1l  ^  mass  transportation  effect. 

This  average  elevation  of  the  water  level  can  be  roughly  estimated 
by  computing  the  discharge  due  to  a  succession  of  limit  solitary 
waves  at  the  entrance.  Then,  by  assuming  that  this  volume  of  water 
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is  spread  northward  In  the  bay  at  a  speed  Vgd”,  this  permits 
calculation  of  the  total  horizontal  area  on  which  an  elevation  of 
the  mid-water  level  exists  by  mass  transportation.  It  is  known 
that  the  mass  transportation  due  to  a  solitary  wave  is; 


II 


{ft. -sec.  units) 


where  d  is  the  depth  and  H  is  the  height.  Since  in  the  case  of 
a  limit  solitary  wave  H  =  0.78  d,  the  total  volume  of  water  enter¬ 
ing  the  Chesapeake  Bay  over  a  length  is: 

V  =  1.7N  /  d®  d/ 


where  N  is  the  number  of  waves. 

6 .  Resonance  of  Harbor  Basins 

The  penetration  of  comparatively  small  waves  into  local  har¬ 
bor  basins  could  cause  severe  agitation  if  the  wave  period  is  close 
to  a  resonant  frequency  of  the  harbor  basin.  The  majority  of 
harbor  basins  in  the  Newport/Norfolk  area  consist  of  two  parallel 
piers  open  to  the  main  channel.  The  natural  resonant  frequencies 
in  the  longitudinal  direction  are  phenomena  very  similar  to  the 
r  esonance  of  sound  in  organ  pipes  (nee  Defant  (1961)>  I_»e  Mehaute 
{195d,  I96l,  1962),  Miles  and  Muhk{196l)  )-Thc  amplitude  of 
agitation  in  the  harbor  basin  at  resonance  theoretically  tends  to 
infinity  if  friction  is  neglected.  In  practice  also,  the  agitation  is 
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limited  physically  by  a  maximum  wave  steepness.  In  some  basins 
it  is  also  possible  that  resonance  could  occur  in  a  lateral  direction. 

A  standing  wave  of  length  equal  to  twice  the  width  of  the  basin 
would  give  such  conditions. 

The  resonant  wave  frequency  in  the  lateral  direction  is  given 
by  T  =  where  b  is  the  width  of  the  basin  and  ,y/gd  Is  the 

velocity  of  propagation  of  a  long  wave  in  shallow  water.  In  fact 
there  could  also  exist  harmonics  but  it  will  be  seen  that  only  the 
fundamental  lateral  resonant  frequency  is  of  interest  here. 

The  longitudinal  resonance  phenomenon  is  a  little  more  com¬ 
plicated.  When  the  harbor  basin  is  long  and  narrow,  the  resonant 

4  ^ 

wave  lengths  are  given  by  L^  s  ^  j  where  V  is  the  length  of 

the  basin  and  n  is  an  integer  (0,  1,  Z,  3.  .  L  =  T^  gd  gives 
4  ^ 

(Xn  -t-  1 )  resonant  wave  periods.  In  practice  when 

the  width  of  the  harbor  basin  is  not  negligible  with  respect  to  the 
length  of  the  basin,  a  correction  factor  has  to  be  applied  such  that: 

Tn  -  T„(l  +  e) 

where  1  +  C  Is  a  function  of  b/^  and  L/V  (i.e.  with  L/-/  =  0  ). 
The  resonance  phenomenon  can  be  closely  evaluated  from: 


where  X  is  Euler's  constant  t  0.577215).  This  equation  has  to 
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be  solved  by  trial  and  error  for  0  .when  b/-/  is  given  for 
various  values  of  n.  Two  graphs  which  greatly  facilitate  this 
procedure  are  presented  as  Figs.  5  and  6.  The  general  method 
of  obtaining  the  resonant  periods  can  be  described  as  follows. 

(1)  For  a  given  value  of  n  corresponding  to  the  harmonic  to  be 
investigated  a  value  of  0  has  to  be  assumed.  This  value  can 


be  estimated  quite  closely  with  practice.  (E)  Fig.  5  is  used  to 

/3  b  _ 

derive  /  „  o\  (3)  /  can  be  computed  from 

H’-f)  '('-¥■) 

the  basin  characteristics.  (4)  Using  this  value  of  /  „  q\ 


'('-¥) 


the  value  of  the  right  hand  term  above  can  be  read  off  Fig.  6. 

P 

(5)  This  value  should  check  with  - - ~:x'  calculated  ii 

,  4  1-f.) 

Step  Z.  \  ^  / 


calculated  in 


In  practice  it  is  found  that  the  correction  term  f  is 
almost  negligible  for  any  resonant  freqtiency  but  the  fundamental 
(n  J=  0).  For  harmonics  (n  ^  D),  0  may  be  considered  as  a 

constant  within  the  usual  range  of  as  is  seen  in  Fig.  6, 

The  amplitude  of  agitation  at  resonance  has  been  estimated 
from  an  adaptation  of  a  theory  and  some  experimental  results  by 
Le  Mehaute  (1959,  1961)  for  the  agitation  in  a  harbor  caused  by 
waves  at  least  four  or  more  times  the  wave  amplitude  in  the  main 
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/3  =  L/  I 

FIGURE  5 

NOMOGRAPH  FOR  HARBOR  RESONANCE 
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FIGURE  6 

NOMOGRAPH  FOR  HARBOR  RESONANCE 
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CHAPTER  V 

MAIN  QUANTITATIVE  RESULTS 

The  rnoat  important  results  have  already  been  summarized  in 
Chapter  II,  Part  A  of  this  report.  In  particular.  Fig.  3  representing 
a  map  of  the  maximum  possible  flooded  areas,  gives  the  essential  re¬ 
sults  of  this  study. 

The  purpose  of  this  chapter  is  to  complement  the  results  already 
given  by  presenting  more  information  and  results  given  by  the  interme¬ 
diate  steps  of  the  calculations  without  entering  into  all  the  details  pre¬ 
sented  in  the  various  appendices .  Those  quantitative  results  have  been 
obtained  by  applying  the  methods  of  calculation  described  in  Chapter  IV 
to  the  Atlantic  coast,  with  special  considerations  for  the  Chesapeake  Bay 
area. 

The  original  disturbance  was  assumed  to  be  located  at  a  point 
defined  by: 

Longitude  35^  57-1/2'  N 

Latitude  68°  25'  W 

where  the  depth  is  3,000  fathoms,  but  surrounded  by  a  depth  of  2,600 
fathoms.  (Fig.  2) 

The  wave  characteristics  were  provided  for  three  values  of  the 
magnitude  of  the  disturbance  at  20  miies  from  their  origin.  They  are 
defined  by  using  an  asymptotic  solution  developed  by  Kranzer  and  Keiler, 
which  is  assumed  to  be  valid  at  a  sufficient  distance  from  the  disturbance 
(See  Appendices  I  and  11).  The  wave  input  is  characterized  by  wave  trains, 
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the  distance  between  crests  varying  with  time  as  shown  by  Fig,  1. 

Considering  the  particular  location  o£  the  disturbance  and  the  bottom 
contours,  it  appears  that  refracti-^n  effects  are  absolutely  negligible  except 
at  some  locations  near  the  shoreline  and,  in  particular,  at  the  entrance  to 
Chesapeake  Bay,  The  bottom  contours  are  practically  concentric  circles 
around  the  location  of  the  disturbance  as  shown  by  Fig,  7.  Hence,  up  to 
the  continental  slope,  the  plan  of  the  wave  crests  is  represented  by  a  series 
of  concentric  circles,  despite  any  shoaling  effect.  The  refraction  effects 
over  the  continental  shelf  appear  only  locally  and  can  be  neglected.  It  is 
for  this  reason  that  the  distance  between  ortbogonals  given  in  Appendix  V 
can  be  replaced  by  a  linear  function  of  the  distance  from  the  origin  R. 

It  is  for  this  reason,  also,  that  a  wave  refraction  diagram  has  not  been 
included  since  it  would  simply  give  the  radii  as  orthogonals  even  for  the 
first  longest  waves  which  ''feel”  the  bottom  at  an  early  stage. 

The  results  of  all  calculations  are  presented  in  the  enclosed  tables 
and  graphs.  Fig.  7  shows  the  bottom  profile  from  the  original  disturbance 
up  to  the  continental  slope  along  four  significant  radii  shown  on  Fig.  2. 

This  figure  demonstrates  the  relative  degree  of  accuracy  obtained  by 
assuming  that  the  bottom  contours  are  concentric  circles. 

Fig.  8  summarizes  the  main  results  which  have  been  obtained  by 
applying  the  methods  of  calculation  presented  in  Appendix  II  for  the  deep 
water  wave  problems,  A  correction  factor  for  the  first  longest  wave  has 
also  been  calculated  to  take  into  account  the  shoaling  effect.  The  results 
have  been  calculated  for  a  radius  of  300  nautical  miles,  i.e.  the  toe  of 
the  continental  slope , 
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FIGURE  8 

WAVE  PERIOD  AND  WAVE  ENVELOPES  VERSUS  TIMES  AT 
300  NAUTICAL  MILES  (TOE  OF  THE  CONTINENTAL  SLOPE- 
DEPTH  6000  FEET) 
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Thti  coefficient  of  transmission  over  the  continental  slope  has 
been  calcuiatei)  as  shown  in  Chapter  IV.  The  values  of  this  coefficient 
of  transmission  depend  upon  the  choice  of  depths  d^  and  d^  and  the 
wave  period.  The  period  is  well  known  and  varies  with  time  as  shown 
in  the  previous  graphs  (Fig,  8).  But  the  choice  of  depths  is  quite  arbi¬ 
trary  because  the  bottom  does  not  present  any  sharp  angle  which  would 
permit  a  clear  interpretation  of  d^  and 

Considering  the  bottom  profile  presented  in  Fig.  7,  it  seems 
that  the  most  reasonable  value  for  is  1000  fathoms,  which  for  most 
of  these  wave  periods  can  be  considered  mathematically  as  infinity. 

The  continental  slope  ends  at  a  depth  of  d^  =  ZO  fathoms,  but 
in  practice  the  slope  between  120  fathoms  and  20  fathoms  is  very  gentle. 
Hence  it  can  be  assumed  that  all  the  wave  energy  is  transmitted  from  a 
depth  of  120  fathoms  to  a  depth  of  20  fathoms,  and  the  wave  reflection 
ia  due  to  that  part  of  the  slope  below  a  depth  of  120  fathoms  only. 

However,  in  order  to  be  able  to  estimate  the  error  which  is  intro¬ 
duced  by  choosing  other  combinations  of  depths  (and  more  advanced  theory), 
the  coefficient  of  transmission  has  also  been  calculated  for  two  extreme 
cases,  namelyr  (1)  all  wave  energy  is  transmitted  over  thcj  continental 
shelf  (there  is  no  reflection),  and  (2)  the  continental  shelf  acts  as  a 
vertical  wall  from  dj  ~  1,000  fathoms  to  =  20  fathoms.  The  values 
obtained  in  the  intermediate  case  considered  as  the  best  case,  and  these 
two  extreme  cases,  arc  presented  in  Tabic  1  as  a  function  of  the  wave 
period. 

It  is  seen  that  except  for  the  first  longest  waves  (which  are  not  the 
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TABLE  1 

COEFFICIENT  OF  TRANSMISSION  OVER  THE  CONTINENTAL  SHELF 
VERSUS  WAVE  PERIOD 


Wave  Period 

T  =  sec. 

No  Reflection 

Reflection  from 
1000  fath.  to 

120  fath. 

T  ransmission 
from  120  fath. 
to  20  fath. 

Reflection 
from  1000  fath. 
to  20  fath. 

343 

3.  75 

2.55 

1.98 

257 

3.25 

2.  40 

1.95 

188 

2.  80 

2.  31 

1.87 

171 

2.  65 

2.25 

1.75 

147 

2.  45 

2.  16 

1.72 

128 

2.  31 

2.09 

1.69 

114 

2.  17 

2.02 

1.65 

103 

2,07 

1.96 

1.62 

85.8 

1. 88 

1.83 

1.56 

73.5 

1. 74 

1.71 

1.60 

64.4 

1.64 

1.61 

1.46 

57.2 

1.55 

1.53 

1.41 

51.5 

1. 47 

1.46 

1.37 

46.8 

1.40 

1.39 

1.32 

42.9 

1. 37 

1.36 

1.29 

39.  6 

1 . 35 

1 .  35 

1.26 

38.8 

1.21 

1.  31 

1.23 

34.3 

1.27 

1.27 

1.20 

32.2 

1.22 

1. 22 

1.17 

30.3 

1.19 

1.19 

1 .14 

Z8 . 6 

1.16 

1.16 

1 . 10 

27.1 

1.11 

1.11 

1.09 

25.7 

1.09 

1.09 

1.07 

24.  5 

1.07 

1.06 

1.06 

23.4 

1 , 05 

1.05 

1.05 
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highest  waves)  the  difference  is  very  small  and  consequently  the  error 
due  to  the  proposed  approximations  is  quite  negligible.  It  is  also  seen 
that  the  waves  of  the  second  wave  train  are  relatively  so  short  that  they 
pass  over  the  shelf  without  being  significantly  ii\fluenced  by  the  conti¬ 
nental  slope,  (A  more  exact  theory  will  be  necessary  for  tsunami  waves 
due  to  seismic  disturbances  because  being  of  longer  wave  period,  these 
waves  are  more  influenced  by  the  continental  slope.) 

Now  the  wave  trains  over  the  continental  shelf  can  be  calculated 
by  a  simple  multiplication  of  the  coefficient  of  transmission  with  the  wave 
heights  given  at  ^00  miles  by  taking  into  account  the  variation  of  wave 
period  with  time  and  the  spreading  effect.  Hence  the  wave  trains  at  320 
miles  from  the  explosion  over  a  depth  of  20  fathoms  are  given  in  Fig,  9, 

The  wave  motion  over  the  continental  shelf  no  longer  depends  upon 
the  period  and  shallow  water  wave  theories  can  be  applied.  The  bottom 
friction  is  of  great  importance.  Then  the  calculation  has  to  be  continued 
step  by  step  as  indicated  in  Chapter  IV.  The  results  of  such  calculations 
are  presented  in  Table  2  and  in  Fig.  10.  The  interval  AR  has  been 
chosen  to  be  five  miles.  The  maximum  wave  heights  corresponding  to 
the  three  cases  under  study  are  also  indicated  on  the  graphs, 

It  is  seen  that  in  Case  II  and  Case  TII  most  of  the  waves  reach  the 
breaking  index  curve  prior  to  reaching  a  bottom  slope  smaller  than  1/100 
where  the  breakers  become  saturated.  Hence  the  waves  follow  the  breaking 
index  curve  whatever  their  original  wave  height  in  deep  water.  It  is  seen 
that  only  in  Case  I  does  the  wave  reach  the  breaking  index  curve  over  a 
slope  which  is  steep  enough  for  the  breaker  to  be  saturated.  Hence  there 
is  less  wave  run-up  than  in  Cas^s  11  and  III, 
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FIGURE  9 

WAVE  PERIOD  AND  WAVE  ENVELOPES  VERSUS  TIMES  AT 
320  NAUTICAL  MILES  (EDGE  OF  CONTINENTAL  SHELF) 
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A  very  simple  conclusion  can  be  deduced  from  these  calculations: 
the  maximum  possible  wave  run-up  is  the  same  in  Cases  II  and  III  for 
practically  all  the  waves.  This  maximum  wave  run-up  is  about  twice  the 
value  of  the  breaking  wave  height  where  the  slope  is  O.Ol.  Such  a  slop 
is  encountered  on  the  lower  part  of  the  entrance  of  Chesapeake  Bay  over 
a  depth  varying  from  2.0  feet  to  30  feet.  Hence  the  maximum  vertical 
ruU“Up  over  an  average  slope  of  I /I S  is  ^  x  0 . 78  x  ft.  =  feet. 

It  can  therefore  be  concluded  that  the  dunes  which  are  only  15  to  20  feet 
above  the  sea  level  will  certainly  be  overtopped  and  breached.  These  dunes 
once  overtopped,  are  very  rapidly  destroyed  by  ei'osion.  This  destruction 
occurs  with  the  first  few  waves.  Then  the  following  waves  penetrate 
deeper  and  deeper  inshore  destroying  all  existing  coastline  structures  by 
dynamic  action. 

However,  it  can  be  expected  that  this  dynamic  effect  will  be  limited 
to  the  neighborhood  of  the  shoreline  because  the  waves  are  rapidly  damped 
as  they  penetrate  further  over  land.  Since  the  coastline  is  very  flat,  a 
large  area  can  be  flooded  by  a  rise  in  water  level  caused  by  the  mass  trans¬ 
portation  effect  due  to  the  series  of  "solitary"  waves. 

On  the  northern  part  of  the  entrance  of  Chesapeake  Bay  the  bottom 
slope  has  a  small  value  before  reaching  the  shoreline.  Hence  the  natural 
protection  of  the  shelf  is  extended  almost  up  to  the  shoreline.  The  waves 
lose  their  energy  over  the  existing  swamps  and  give  a  relatively  smaller 
run-up  on  dry  land. 

At  the  entrance  of  the  Bay  itself  the  wave  heights  in  Cases  II  and 
in  are  directly  related  to  the  depth  by  the  relationship  H  =  0. 78  d.  This 
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very  simple  law  results  from  the  fact  that  the  waves  follow  the  breaking 
index  curve  and  appear  at  the  entrance  of  Chesapeake  Bay  as  a  succession 
of  almost  symmetrical  spilling  breakers  which  can  be  considered  as  a 
succession  of  limit  solitary  waves.  Since  the  average  depth  is  ?.6  feet 
(including  2, 8  feet  of  tide),  they  will  appear  as  a  succession  of  solitary 
waves  of  about  20  feet  in  height. 

Further  quantitative  information  has  been  obtained  in  Appendix  VI 
for  a  i/lQ  bottom  slope.  Despite  the  fact  that  the  slope  is  usually  smaller, 
these  values  give  a  rough  estimate  of  what  can  be  expected  in  the  present 
case , 

First,  it  has  been  seen  that  the  height  of  the  leading  front  when 
reaching  the  shoreline  on  the  open  coast  is  about  seven  feet  (  -7  /d^  =  0.2E 

and  d^  =  30  feet).  The  thickness  of  the  water  which  follows  the  leading 
front  will  increase  up  to  a  value  of  about  13  feet  {see  Fig.  VI-21).  The 
speed  of  the  water  of  the  leading  front  at  the  shoreline  is  48.6  ft/sec 
(  Ug  /  '/^^i  ~  1'568). 

The  maximum  horizontal  run-up  for  the  leading  front  (assuming 
this  leading  front  is  not  caught  up  by  the  following  wave  elements)  gives  the 
minimum  horizontal  distance  where  the  waves  may  have  a  dynamic  effect. 
This  value  found  by  giving  a  zero  value  to  u  in  the  formula 

(VI-56).  Then  it  is  found  that  could  be  as  high  as  1000  feet. 

In  fact  the  bottom  slope  is  so  small  after  the  dune  has  been  over- 
topped  that  the  obtained  formula  gives  only  the  very  minimum.  Because 
of  the  limitation  of  validity  of  the  obtained  formula,  it  must  be  expected 
that  the  dynamic  effect  of  the  succeeding  wave  will  be  spread  over  a  3000 
foot  area  from  the  shoreline.  Also,  since  the  water  depth  will  increase 
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because  of  the  mass  transportation  over  the  dune,  the  wave  effect  will  bo 
felt  deeper  and  deeper » 

The  present  state  of  the  art  is  such  that  it  is  not  possible  to  give 
presently  more  information  about  this  problem,  which  would  require  taking 
account  of  the  dune  evolution  by  erosion.  Fig.  3  gives  what  can  be  con- 
sidered  as  the  maximum  possible  flooding  areas  which  will  certainly  not 
be  exceeded. 

Fig.  n  presents  the  wave  refraction  diagram  for  the  spilling 
breakers  entering  the  Chesapeake  J3ay.  It  is  seen  that  the  wave  motion 
rapidly  becomes  very  complex.  A  number  of  orthogonal?  cross  each 
other  limiting  the  value  of  the  method,  in  spite  of  the  improvements  for 
nonlinear  effects  proposed  in  Chapter  IV.  It  is  interesting  to  note  the 
differences  in  wave  refraction  between  a  spilling  breaker  as  presented  in 
Fig.  11  and  a  bore  as  presented  in  Appendix  VII. 

Also,  it  is  known  that  when  a  long  wave  passes  over  a  shoal,  a 
double  crest  appears.  This  fact  has  been  observed  in  practice  at  sea,  as 
well  as  in  a  laboratory.  The  long  waves  are  transformed  into  a  succession 
of  shorter  waves.  Such  a  phenomenon  can  be  analyzed  only  in  the  case  of 
a  simply  shaped  two-dimensional  shoal.  It  is  impossible  to  analyze  it  for 
a  two-dimensional  case  with  the  present  state  of  the  art. 

Whatever  the  probable  quantitative  accuracy  of  the  motion  calcu¬ 
lated  in  Chesapeake  Bay,  the  wave  run-up  has  been  evaluated  from  the 
wave  refraction  diagram  presented  in  Fig.  11  and  the  methods  presented 
in  Chapter  IV.  The  results  of  these  computations  are  also  given  by  Fig.  3. 

’^Because  of  its  size.  Fig.  J 1  is  included  At  the  back  of  this  report 
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Because  o£  the  complex  wave  patterns  just  inside  the  entrance  to 
Chesapeake  Bay,  it  is  difficult  to  assess  the  damage  to  the  Chesapeake 
Bay  bridge “tvinnel  which  is  currently  under  construction.  Although  the 
mean,  wave  height  is  about  iiO  feet,  because  of  the  formation  of  short - 
crested  waves,  the  local  wave  heights  in  parts  of  the  structure  could 
greatly  exceed  this  value  and  he  as  high  as  40  feet  locally.  Also  there 
could  be  serious  erosion  effects  around  the  foundations.  Knowledge  of 
tho  movement  of  sand  and  silt  by  long  waves  is  very  Umited.  These  same 
considerations  apply  to  the  Newport-Norfolk  bridge -tunnel  from  Willoughby 
Spit  to  Old  Point  Comfort. 

Even  though  the  structures  may  be  saved  from  destruction  by 
dynamic  effects  or  undercutting  of  their  foundations,  any  wave  over¬ 
topping  near  the  tunnel  entrances  would  flood  out  the  tunnel  sections  and 
render  the  structures  inoperative  for  traffic. 

The  mid-elevation  in  the  lower  Chesapeake  Bay  due  to  the  mass 
transportation  by  the  waves  at  the  entrance  has  been  found  to  be  about 
one  foot.  In  any  case,  it  cannot  exceed  tw'o  feet.  However,  this  average 
elevation  generates  a  very  long  wave  which  propagates  northward  at  a 
velocity  -v/gd  similar  to  a  tidal  wave.  Hence  the  elevation  in  the  upper 
Chesapeake  Bay  and  Potomac  River  can  be  as  high.  For  example,  the 
tide  is  slightly  higher  at  Washington  than  at  the  entrance  to  Chesapeake 
Bay.  However,  the  shorter  waves  due  to  explosion  are  quickly  damped 
by  successive  diffraction  and  do  not  penetrate  so  far. 

Now  the  results  of  the  problem  of  harbor  resonance  are  presented. 

About  one  hour  after  the  v/aves  arrive  at  the  entrance  to  Chesapeake 
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Bay,  they  will  arrive  at  Newport  News  and  a  short  while  later  their  effects 
will  be  apparent  in  the  Elizabeth  River  at  Norfolk  and  Portsmouth.  The 
fundamental  natural  period  of  oscillation  in  the  longitudinal  direction,  the 
first  two  harmonics,  and  the  natural  period  of  lateral  oscillation  for  various 
harbor  basins  have  been  computed  and  are  presented  in  Tables  3  through 
6.  The  locations  of  these  basins  are  shown  on  Fig.  12  through  14.  It  is 
seen  that  the  range  of  wave  periods  at  which  resonance  will  occur  compares 
very  closely  with  the  range  of  wave  periods  in  the  wave  trains  generated  by 
the  disturbance  as  they  arc  shown  by  Fig.  15. 

Also,  in  some  of  the  harbor  basins  and  navigation  channels  serious 
siltation  and/or  erosion  could  occur  very  rapidly  because  of  the  high 
velocities  and  severe  agitation.  In  particular,  the  navigation  channel  of 
Little  Creek  would  certainly  be  filled  by  sand.  This  will  be  especially  true 
if  resonance-  occurs  as  the  motion  is  then  amplified.  Undermining  by 
scouring  would  occur  at  pier  foimdations  and  the  material  deposited  in 
relatively  calm  water  in  the  main  shipping  channels.  Such  phenomena  have 
already  been  observed  by  the  action  of  tsunami  waves.  (Iwasaki  and 
Horikawa  (I960)  ), 
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TABLE  3 

NATURAL  OSCILLATORY  PERIODS  OF  HARBOR  BASINS  IN  THE 
NEWPORT  NEWS  AREA 

Harbor  Basin 

F  undamental 
Period  of 

First 

Se  cond 

Hate  ral 
Oscillation 

No, 

Oscillation 

Harmonic 

Harmonic 

Pe  riod 

Sec . 

Sec . 

Sec . 

Sec . 

1 

816 

261 

142 

32 

2 

125 

34 

19 

46 

3 

123 

38 

22 

46 

4 

155 

45 

25 

18 

5 

151 

40 

23 

49 

6 

104 

29 

17 

25 

7 

95 

28 

16 

9 

8 

97 

29 

16 

9 

9 

96 

29 

16 

9 

10 

86 

23 

13 

17 

1 1 

159 

46 

26 

30 

12 

146 

42 

24 

3U 

I  3 

153 

44 

25 

30 

14 

129 

38 

22 

20 

15 

129 

38 

22 

20 

16 

101 

27 

15 

19 

17 

117 

32 

18 

16 

18 

96 

27 

15 

13 

ApproxirnHtf’  prrivni  time  after  enterint^  Chesapeake  Bay; 
Slightly  greater  than  one  hour. 


-7^- 

CONFIDENTIAT. 


CONriOf^NTIAL 


TABLE  4 

NATURAL  OSCILLATORY  PERIODS  OF  HARBOR  BASINS 
IN  NORFOLK  HARBOR  REACH 


Harbor  Basin 
No. 

Fundamental 
Period  of 
Oscillation 

First 

t  T_  _  „  1 

1  ifi  1  iiionLc 

Second 

Harmonic 

Lateral 

Oscillation 

Period 

Sec . 

Sec . 

Se  c . 

Sec . 

19 

321 

96 

55 

1  19 

20 

400 

102 

58 

190 

21 

203 

50 

34 

21 

22 

191 

54 

31 

24 

23 

19i 

54 

31 

24 

24 

196 

58 

33 

19 

25 

196 

58 

33 

19 

26 

156 

45 

26 

18 

27 

379 

105 

60 

60 

28 

122 

33 

19 

27 

29 

121 

32 

18 

32 

30 

173 

47 

27 

33 

31 

378 

113 

64 

32 

32 

230 

59 

34 

78 

33 

261 

74 

42 

35 

34 

188 

54 

32 

24 

Approximate  arrival  time  after  entering  Chesapeake  Bay; 


Slightly  less  than  one  hour. 
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FIGURE  13 

HARBOR  BASINS  IN  NORFOLK  HARBOR  REACH 
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TABLE  5 


NATURAL  OSCILLATORY  PERIODS  OF  HARBOR  BASINS  IN  PORT 
NORFOLK  AND  TOWN  POINT  REACH  ON  THE  EAST  SIDE  OF 
ELIZABETH  RIVER 


Harbor  Basin 
No. 

F  undamental 
Period  of 
Oscillation 

first 

Harmonic 

Second 

Harmonic 

Late  ral 
Oscillation 
Period 

Sec. 

Sec. 

Sec , 

Sec. 

35 

165 

46 

26 

24 

36 

191 

54 

31 

24 

37 

167 

47 

27 

23 

38 

172 

50 

29 

20 

39 

165 

48 

27 

18 

40 

186 

54 

31 

19 

41 

183 

54 

31 

16 

42 

68 

19 

11 

9 

43 

68 

19 

11 

9 

44 

68 

19 

11 

9 

45 

104 

32 

18 

7 

46 

113 

33 

19 

12 

47 

159 

50 

28 

8 

48 

no 

33 

19 

8 

49 

93 

28 

16 

8 

50 

92 

26 

15 

14 

51 

126 

36 

20 

15 

52 

65 

18 

10 

10 

53 

113 

yi 

19 

12 

Approximate  arrival  time  after  catering  Chesapeake  Ray: 
1-1/4  hours , 
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TABLE  6 

NATURAL  OSCILLATORY  PERIODS  OF  HARBOR  BASINS  IN  THE 
LOWER  REACH  AND  WEST  SIDE  OF  PORT  NORFOLK  REACH  IN 
ELIZABETH  RIVER 


First 

Ha  rmonic 

Second 

Harmonic 

Lateral 

Oscillation 

Period 

Sec. 

Sec. 

Sec. 

Sec. 

54 

59 

17 

10 

8 

55 

25 

14 

14 

56 

115 

32 

18 

18 

57 

114 

33 

18 

13 

58 

102 

29 

17 

12 

59 

112 

32 

18 

13 

60 

121 

32 

19 

27 

61 

1 18 

32 

18 

26 

62 

118 

32 

18 

26 

63 

118 

32 

18 

26  . 

64 

142 

40 

23 

17 

65 

142 

40 

23 

17 

66 

142 

40 

23 

17 

67 

62 

18 

10 

7 

68 

129 

40 

22 

10 

69 

72 

21 

12 

8 

70 

79 

21 

12 

15 

71 

88 

26 

15 

9 

72 

104 

31 

17 

9 

73 

142 

44 

25 

8 

74 

142 

44 

25 

8 

75 

142 

44 

25 

8 

76 

i  42 

44 

4b 

B 

77 

142 

44 

25 

12 

78 

267 

70 

40 

67 

Approximate  arrival  time  after  entering  Chesapeake  Bay: 
1-1/4  hours . 
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0  5  0  100  160  200  250  300  350 

TIME  IN  MINUTES 


FIGURE  15 

WAVE  PERIOD  VERSUS  TIME  AT  THE  ENTRANCE  OF 
HARBOR  BASINS 
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CHAPTER  VI 

ON  THE  ACCURACY  OF  THE  OBTAINED  RESULTS  AND 
RECOMMENDATIONS  FOR  FURTHER  INVESTIGATIONS 

1 .  The  theory  of  Kranzer  and  Keller  and  its  application  to  the  pro¬ 
blem  of  the  distiirbance  caused  by  a  nuclear  explosion  is  probably  the 
best  approach  presently  available  to  analyze  the  phenomenon.  However, 
because  of  the  simplifying  assumptions  necessary  to  relate  the  factvial 
disturbance  to  the  mathematical  model  some  inaccuracy  will  be  inherent 
in  the  input  data  as  shown  in  Appendix  I.  It  is  believed  that  the  presented 
results  of  this  report,  whatever  the  difficulty  encountered  in  their  calcu¬ 
lation,  are  no  less  accurate  --  and  probably  more  accurate,  than  the 
input  definition.  It  must  be  understood  that  this  statement  is  not  a  crit¬ 
icism  of  the  application  of  the  Kranzer  and  Keller  theory  for  defining  the 
input  because  no  better  method  presently  exists. 

Hence,  if  further  investigations  are  considered,  the  particular 
problem  of  wave  agitation  dvie  to  explosion  should  be  one  of  the  problems 
of  primary  importance.  This  particular  problem  can  be  subjected  to 
field  experimentation.  It  can  also  he  handled  by  further  theoretical  anal¬ 
ysis,  First  a  mathematical  model  for  the  impulse  needs  to  be  related  to 
the  yield  of  the  nuclear  explosion.  Then  the  wave  motion  due  to  this  im¬ 
pulse  definition  (which  can  be  different  from  that  chosen  by  Kranzer  and 
Keller)  can  be  calculated. 

2.  The  Kranzer  and  Keller  theory  has  been  established  for  a  hori¬ 
zontal  bottom.  The  effect  of  a  change  of  depth  is  also  an  important  aspect 
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which  should  require  further  analysis  with  the  aim  of  determining  the 
effect  of  the  continental  slope.  This  problem  can.  be  solved  by  modifi¬ 
cation  of  the  Kranzor  and  Keller  theory.  In  the  case  of  a  bottom  slope, 
the  Bessel  fvinctions  should  be  replaced  by  the  Mathieu  functions.  The 
problem  can  also  be  approached  by  analyzing  the  problem  of  long  waves 
over  a  bottom  slope.  This  work  then  should  begin  with  a  literature  sur¬ 
vey  and  synthesis  of  the  existing  theories  on  this  subject.  This  problem 
can  also  be  studied  quite  easily  experimentally  in  a  two-dimensional 
wave  tank^  specially  adapted  for  the  problem. 

3,  Now  considering  the  problem  of  the  wave  deformation  over  the 
continental  shelf,  it  can  be  said  that  tlwre  is  presently  no  satisfactory 
theory  which  takes  into  account  all  the  important  factors  such  as  the 
bottom  friction  and  the  vertical  acceleration  component  (see  Morrison 
and  Crook  (1954)  ).  There  is  a  need  for  a  the<iry  for  long  waves  over  a 
gentle  slope  where  the  vertical  acceleration  tern>  and  path  curvature  ef¬ 
fects  are  taken  into  account.  Then  the  determination  of  the  friction  co¬ 
efficient  must  rely  on  carefxil  field  experiments.  If  one  neglects  the 
bottom  friction  term,  a  possible  approach  to  the  problem  exists  in  neg¬ 
lecting  the  quadratic  term  w  — —  in  the  Euler  equation,  which 

gives  in  the  usual  symbolism: 

b  u  b  u  .  b  p 

bt  ^  bx"  p  h  X. 

along  a  horizontal  axis  OX.  AU  quadratic  terms  along  the  vertical  axis 
OZ  could  be  neglected,  namely  u  — ;  —  ,  w  — ^  -  —  ,  such  that 

o  X  o  a 

b  w  _  1  b  (p  P  g  z) 

8  t  p  h  z 
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In  a  word,  the  motion  is  nonlinear  horizontally,  but  iinearlaed  vertically. 
However,  it  is  pointed  out  that  in  the  usual  theory  of  long  waves, 

o  t 

is  neglected  and  the  pressure  is  assumed  to  be  hydrostatic.  Such  an 
assumption  gives  a  wave  deformation  even  on  a  horizontal  bottom.  This 
deformation  is  never  observed  in  practice  on  a  gentle  slope  such  as  that 
of  the  continental  shelf.  The  neglect  of  ^  leads  to  the  Earnshaw 

paradox  discussed  in  Appendix  VI. 

Combining  these  two  equations  with  the  continuity  and  boundary 
conditions  (Wj^  =  S  at  the  bottom  of  slope  S  ),  a  solution  can  be  found 
by  assuming: 


a.  w  is  Linearly  distributed  from  the  bottom  (at  depth  d) 
to  the  free  surface  (at  elevation  ^  above  the  mean  water  level): 


.  d  +  z 

w  =  w,  4  —r-c — ^ 
b  <1  +  17 


-42L_ 

b  t 


b.  In  first  approximation  u  can  be  assumed  to  be  a  constant 
over  a  vertical:  u  -  u  (x,t)  .  Then,  to  second  and  third  approximation: 

u  -  u  (x,t)  I  +  A  (d  4  k)  4  B  (d  + 

The  problem  then  consists  of  determining  the  coefficients  A  and  B.  A 
bottom  friction  term  can  also  be  added,  but  a  good  determination  of  the 
friction  factor  would  require  analysis  of  field  results. 


4.  A  theoretical  approach  to  wave  run-up  has  been  developed  in  this 

report  (Appendix  VI}.  It  can  be  said  that  the  run-up  can  now  be  calculated 
theoretically.  However,  in  order  to  obtain  practical  results  from  the 
theory,  further  studies  remain  to  be  done,  namely  the  calculation  of  the 
wave  run-up  for  various  slopes  and  various  friction  coefficients.  This 
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can  be  done  after  development  of  a  computing  program  which  would  be 
xiecessary  because  of  the  length  of  the  calculations.  Some  experiments 
in  a  specially  built  wave  flurno  for  long  waves  will  also  be  of  great  value 
in  order  to  verify  the  theory  presented  and  give  some  information  on 
the  friction  coefficient.  Analysis  of  field  tests  will  also  yield  a  major 
contribution  which  will  complement  the  theoretical  research. 

5.  A  new  numerical  method  for  studying  the  problem  of  two-dimensional 
waves  in  bays  and  estuaries  has  been  developed  in  this  report  (Appendix 
VII).  It  has  also  been  applied  to  the  lower  Chesapeake  Bay.  However, 

this  method  can  render  more  information  after  an  extension  of  the  cal¬ 
culations  over  a  larger  area.  The  accuracy  can  be  improved  by  using 
smaller  intervals. 

One  theoretical  problem  which  must  also  be  investigated  in  con¬ 
nection  with  the  propagation  of  two-dimensional  waves  in  estuaries  is 
the  problem  of  successive  diffraction  of  nonlinear  long  waves.  Practically 
nothing  is  known  in  hydrodynamics  on  nonlinear  diffraction.  This  study 
can  be  approached  by  the  ’’method  of  wave  derivatives”  as  presented  in 
Appendix  VII,  However,  it  would  be  desirable  to  improve  the  presented 
method  by  taking  into  account  the  effects  due  to  vertical  acceleration  and 
and  path  curvatures  as  it  is  prepared  in  Appendix  VI.  Otherwise,  an 
unrealistic  wave  profile  deformation  will  b(3  obtained  in  view  of  the  Karnshaw 
paradox. 

6.  The  problem  of  resonance  in  a  harbor  basin  due  to  periodic  waves 
has  already  been  subjected  to  extensive  theoretical  studies  (McNown  (1951)  ), 
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Le  Mehaut^ ( 1 954,  i96i)»  Miles  and  Munk  (19<iil)  )  .  However,  the  waves 
due  to  explosion  are  not  periodic.  Hence  the  problem  of  harbor  response 
to  a  nonperiodic  wave  must  also  be  analyTied.  It  ’s  a  problem  similar  to 
the  problem  of  harbor  response  to  a  wave  spectrum. 

The  determination  of  the  duration  for  establishing  a  given  state  of 
agitation  under  a  periodic  wave,  or  the  damping  of  this  agitation  when  the 
outside  agitation  ends  ,  are  similar  problems.  They  ail  must  be  solved  by 
equations  of  energy.  Systematic  experiments  in  a  laboratory  will  be  of 
great  help  in  solving  these  problems. 

The  pi'oblom  of  the  response  of  moored  ships  to  seiche  motion  is 
also  one  of  the  topics  which  can  be  considered  for  further  analysis  (Wilson 
( 1951 ,  1959)  ).  Then  it  can  be  reasonably  thought  that  a  solution  of  nmoringa 
could  be  found  which  would  avoid  the  danger  of  wrecking  the  ships  by  col¬ 
lision  against  the  quays  or  against  each  other. 

7.  Finally,  full  use  can  be  made  of  scale  model  technology  and  the 

analog  conxputer.  An  analog  conipviter  is  particularly  interesting  for 
studying  the  propagation  of  linear  periodic  waves.  However,  one  cannot 
expect  too  much  accuracy  by  this  method  because  long  waves  are  nonlinear 
and  the  problem  of  bore  propag-ation  --  which  is  so  important  for  this  pro- 
bleni  --  is  essentially  a  nonlinear  phenomenon,  A  scale  model  automatically 
takes  care  of  nonlinear  effects  and  fa.ithfully  reprodvices  the  bore  phenomenon. 
The  bottom  friction  can  even  be  adjusted  oxperirT«entally. 

Considering  all  the  problems  involved  in  this  study  and  the  size  of 
the  area  to  be  represented,  it  would  be  particularly  convenient  to  have 
several  different  n^odels  ill  various  scales,  A  first  model  at  a  relatively 
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small  horizontal  scale  (say  X  =  i/2»000J  will  represent  the  lower 
Chesapeake  Bay.  A  wave  paddle  could  be.  set  up  at  the  entrance  of 
Chesapeake  Bay  to  reproduce  waves  of  limit  steepness  with  a  wave 
period  which  is  slowly  decreased  with  time  according  to  the  theory  of 
waves  due  to  an  explosion  as  shown  by  Fig.  15.  This  coviid  be  achieved 
by  a  piston  wave  paddle  in  front  of  a  small  gentle  slope.  Such  a  model 
would  need  to  be  distorted  in  order  that  the  depth  would  not  be  smaller 
than  1"  and  averaging  2".  Hence  the  scale  model  would  have  a  vertical 
scale  of  abovit  p.  -  1/200,  i.e.  a  distortion  of  iO.  The  time  scale 

T  "  - -  is  0.007.  Hence  an  HO  second  wave  period  would 

be  represented  by  a  0.48  second  wave,  which  is  above  the  limit  of 
capillarity  wav<is.  Because  of  this  distortion,  this  nmdei  would  have  a 
steeper  slope  than  encountered  in  the  field.  Hence  the  reflection  co¬ 
efficient  may  have  to  be  adjusted  by  use  of  screen  and  roughness  on  the 
concrete  of  the  model. 

The  area  to  be  represented  would  bo  roughly  150,000  x  100,000 
ft.  Efficient  wave  absorbers  would  have  to  be  located  at  the  limit  of 
the  model  because  the  northern  part  of  Chesapeake  Bay  would  not  be  re¬ 
presented.  Hence,  this  indicates  a  75  ft.  x  50  ft.  model. 

It  is  evident  that  all  the  numbers  presented  in  this  discussion  of 
model  stvidies  are  only  for  indicating  the  order  of  magnitude  of  the  model 
size.  Such  a  model  is  possible  even  though  unusual,  but  so  is  the  problem 
under  study.  It  can  reasonably  be  considered  as  the  best  tool  for  solving 
the  problem  within  the  estuary.  A  model  study  as  proposed  would  yield 
the  following  results; 
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a»  The  complete  wave  pattern  in  the  Chesapeake  Bay 
entrance  would  be  faithfully  reproduced,  indicating  bore  phenomena  and 
local  areas  of  high  amplitude,  short ~crested  waves  of  particular  interest 
for  the  Chesapeake  Bay  bridge -tunnel  structure, 

b.  It  is  comparatively  easy  to  make  modifications  on  the 
model  to  determine  the  most  efficient  and  economical  solutions  to  reduce 
damage.  The  best  solutions  can  be  found  by  relative  comparison  of  the 
results,  whatever  the  exactness  of  thair  absolute  value, 

c.  Areas  near  the  shoreline  or  off-shore  which  have  their 
own  special  prohlen^s  will  be  indicated  and  more  detailed  investigations 
can.  be  concentrated  on  those  places. 

d.  Use  of  gilsonite  or  other  movable  material  will  indicate 
areas  where  serious  deposition  of  silt  will  be  likely  to  occur  and  also 
possibly  indicate  areas  subject  to  erosion,  {Le  Mehaut/  and  Collins  (1961)) 

The  probletn  of  resonance  in  a  harbor  requires  a  larger  scale 
model,  but  a  smaller  area  needs  to  be  represented.  Such  a  model  would 
follow  after  the  previous  one  once  the  wave  amplitude  in  the  main  navi¬ 
gation  channels  is  known.  This  model  would  indicate  the  time  of  response 
of  the  harbor  under  consideration  and  the  amplitude  of  the  resonance 
motion. 

These  two  models  would  be  built  after  a  small  preliminary  theore¬ 
tical  vStudy  on  distortion,  dote rm ination  of  the  scales  developed  from  an 
adaptation  of  the  ideas  already  expressed  by  Biesel  and  Le  M^haut^  {1954. 
1955).  Since  many  models  for  long  waves  have  been  built  following  the 
presented  rules:  Alger,  Table  Bay,  Funchal,  and  so  on,  it  seems  that 
they  have  been  successfully  exploited, 
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8,  Finally,  this  problem  cannot  be  handled  with  complete  satisfaction 

until  some  field  tests  arc  performed  to  substantiate  theoretical  research 
and  laboratory  experiments.  Then  a  careful  analysis  of  field  records 
should  give  a  complement  of  information  of  great  value,  f  or  this  purpose 
it  would  be  desirable  to  record  wave  trains  not  only  near  the  shoreline, 
but  also  at  the  edge  of  the  continental  shelf  and  even  offshore  at  the  conti¬ 
nental  slope.  For  this  purpose  a  wave  recorder  needs  to  be  developed 
because  the  present  long  wave  recorders  do  not  operate  satisfactorily  in 
deep  water.  Such  a  problem  is  difficult  and  is  a  challenge  requiring  some 
new  ideas . 
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LIST  OF  SYMBOLS 

a  Horizontal  length  of  the  continental  slope 

A  Function  of  md  (Kq.  2);  also  numerical  coefficient 

b  Width  of  harbor  basin;  also  distance  be.tween  wave 

rays  or  orthogonals 

Distance  between  orthogonals  in  nearshore  zone 

Distance  between  same  orthogonals  at  breaking 
inception 

B  Numerical  coefficient 

C  Wave  (phase)  velocity 

Cy  Value  of  C  in  deepwater  off  the  continental  shelf 

Value  of  C  over  the  continental  shelf 
Chezy  coefficient 
d  Water  depth 

Value  of  d  off  the  continental  shelf;  also  at  any  one 
position;  also  depth  at  which  bottom  slope  is  0.01 

Value  of  d  over  the  continental  shelf;  also  at  any 
second  position 

d^  Value  of  d  at  breaking  inception 

f  Friction  coefficient  (-  g/C^^). 

g  Acceleration  due  to  gravity 

H  Wave  height  from  crest  to  trough 

H  Value  of  H  at  breaking  incention 

b 

1  Length  of  a  harbor  b;^sin 

In  Natural  logarithm  (to  base  e) 

L  Wave  length  between  crests 
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Value  of  L  off  the  continental  shelf 

Value  of  L  over  Ihc  continental  shelf 

Resonant  value  of  L  to  excite  basin  oscillation  of 
n-th  mode 

Wave  number  (=  I  tt/L) 

Value  of  m  off  the  continental  shelf 

Value  of  m  over  the  continental  shelf 

Manning's  friction  coefficient;  also  positive  integer 
(0,  1,  2,  3.  .  .  )  representing  modes  of  oscillation;  also 
numerical  exponent 

Number  of  waves 

Fluid  pressure 

Radial  distance  from  origin  of  disturbance;  also 
run-up  elevation  at  shoreline 

Increment  of  radial  distance 

Bottom  slope  of  sea  bed 

Wave  period 

Period  of  resonance  of  harbor  basin  in  n-th  mode. 
Component  of  water  velocity  in  horizontal  direction 
Value  of  u  at  the  shoreline 

Group  velocity  of  waves  (=  C  for  shallow  water  waves) 
Component  of  water  velocity  in  vertical  direction 
Variable  horizontal  distance  in  a  coordinate  system 
Variable  vertical  distance 

Coefficient  of  transmission  (of  wave  height)  over  the 
continental  slope 

Ratio  {-  L/ y  ) 
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Euler's  constant  (=  0.577215.  .  .  ) 

Shoaling  coefficient 

Elevation  of  free  surface  above  still  water  level 
Horizontal  linear  scale  for  a  model 
Vertical  linear  scale  for  a  rnodel 

Correction  factor  for  computing  oscillation  periods 
Universal  constant  (3.  14159.  .  .  ) 

Mass  density  of  fluid 
Summation  symbol 

Mass  transportation  due  to  a  solitary  wave 
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Defense  Threat  Reduction  Agency 

45045  Aviation  Drive 
Dulles,  VA  20166-7517 


May  6,  1999 


MEMORANDUM  FOR  DEFENSE  TECHNICAL  INFORMATION  CENTER 

ATTN:  OCQ/MR  WILLI7y4  BUSH 


SUBJECT:  DOCUMENT  REVIEW 

The  Defense  Threat  Reduction  Agency's  Security  Offic^ 
has  reviewed  and  declassified  or  assigned  a  new 
distribution  statement: 

-AFSWP-1069,  AD-341090,  STATEMENT  ^  ' 

,^DASA-1151,  AD-227900,  STATEMENT  A  < 

-DASA-1355-1,  AD-336443,  STATEMENT  A^ 

-DASA-1298,  AD-285252,  STATEMENT  A  ^ 


■  DASA-1290, 

-  DASA-1271, 

■  DASA-1279, 

-  DASA-1237, 
•  DASA-1246, 
'  DASA-1245, 

-  DASA-1242, 

-  DASA-1256, 
DASA-1221, 
DASA-1390, 


AD-444208, 

AD-276892, 

AD-281597, 

AD-272653, 

AD-279670, 

AD-419911, 

AD-279671, 

AD-280809, 

AD-243886, 

AD-340311, 


-  DASA- 12^3-,— ad-71  70 97, 


STATEMENT 

STATEMENT 

STATEMENT 

STATEMENT 

STATEMENT 

STATEMENT 

STATEMENT 

STATEMENT 

STATEMENT 

STATEMENT 

STATEMENT 


DASA-1285-5,  AD-443589,  STATEMENT 


A" 

A'^ 

A 
A 
A 
A 
A 
A 
A 

A  — 

A  OK 
A 


'  DASA- 1714, 
■-DASA-2214, 
-■DASA-2627, 
-DASA-2651, 


AD-473132, 
AD-854912, 
AD-514934, 
AD-514615, 
ail-87  669-At- 


STATEMENT 

STATEMENT 

STATEMENT 

STATEMENT 

-&Tj¥TEMENT 


A 

A 

A 

A 

-A- 


DASA-2722T-V3,  AD-518506,  STATEMENT  A 
DNA-3042F,  AD-525631,  STATEMENT  A 
■DNA-2821Z-1,  AD-522555,  STATEMENT  A 


If  you  have  any  questions,  please  call  me  at  703-325- 

1034  . 
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ARDITH  JARRE TT 

Chief,  Technical  Resource  Center 
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